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iteraction of solid state qubits with environmental degrees of freedom strongly affects the qubit 
^~^ dvnamics. and leads to decoherence. In quantum information processing with solid state qubits, 
j^ecoherence significantly limits the performances of such devices. Therefore, it is necessary to fully 
fH understand the mechanisms that lead to decoherence. In this review we discuss how decoherence 
I affects two of the most successful realizations of solid state qubits, namely, spin-qubits and super- 
C^onducting qubits. In the former, the qubit is encoded in the spin 1/2 of the electron, and it is 
^implemented by confining the electron spin in a semiconductor quantum dot. Superconducting de- 

a'ices show quantum behavior at low temperatures, and the qubit is encoded in the two lowest energy 
^evels of a superconducting circuit. The electron spin in a quantum dot has two main decoherence 
■^-^hannels, a (Markovian) phonon-assisted relaxation channel, due to the presence of a spin-orbit in- 
C^eraction, and a (non-Markovian) spin bath constituted by the spins of the nuclei in the quantum 
Cdot that interact with the electron spin via the hyperfine interaction. In a superconducting qubit, 
'yidecoherence takes place as a result of fluctuations in the control parameters, such as bias currents, 
i-^^pplied flux, and bias voltages, and via losses in the dissipative circuit elements. 

o ; 

I Introduction 

^^1.1 What is coherence? and why is it interesting? 

^^'^Coherence is a defining property of quantum mechanics. It can be argued that 
j^quantum coherence is the property that draws a line between the "quantum 
■^world" from the "classical world" . But what exactly is coherence? In physics, 
Q>^tlie term coherence refers to the property of waves to interfere, showing well 
^known interference patterns. Two waves, depending on their relative phase, can 
Q^roduce a constructive interference, characterized by an enhancement of the 
^^mplitude of the wave, or destructive interference, accompanied by a complete 
_ Suppression of it. Only the relative phase of the two waves makes the difference. 
^^To be precise we should therefore speak about phase coherence of quantum 
?-*tates. 

■ - - ' By quantum states here we mean states of a quantum system, which in 
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turn can be constituted by more than one quantum object. The same rules of 
quantum mechanics that allow us to explain and predict interference of one 
object with itself, as the case of an electron through a double slit, predict 
that a system composed by two quantum subsystems can be in a state that 
has no classical counterpart, being a superposition with a precise phase of 
two or more quantum states of the whole system. This property of quantum 
states goes under the name of entanglement, and entangled states need to be 
phase coherent. In particular, coherence, as a property of quantum mechanical 
phenomena, disappears in the classical world, and it is therefore of fundamental 
interest to study it in theory and experiment. It is never completely possible 
to isolate a quantum system from the surrounding world. The system and 
its surrounding interact and, as a result, a randomization of the phase of 
the quantum system takes place, and the initial quantum state ends up in a 
classical state. This process is known as decoherence. 

Only in recent years, thanks to the advances in technology, it has become 
possible to study quantum effects involving single quantum objects, like single 
photons, ions, electron spins, etc... Particular attention has been devoted to see 
coherence, from an experimental point of view, and to understand its limita- 
tion. In fact, though remarkable improvements have been achieved, nowadays, 
to see coherence, a lot of effort must be spent to understand how to preserve 
coherence. In the last decades, the idea of joining quantum physics laws and 
information science gave birth to a new and intriguing branch of science, quan- 
tum information theory, which studies the possibilities that quantum rules offer 
to information processing. In particular, the superposition principle opens the 
possibility to perform new and fast algorithms. The physical implementation 
of quantum information processing represents a challenge because one has to 
deal with the competition between fast and reliable quantum control, that re- 
quires interaction with the outside world, and good isolation of the quantum 
devices in order to ensure long coherence times. 

Therefore, it is important to understand theoretically how decoherence hap- 
pens in the systems under study (here, solid state systems) in order to make 
progress toward this ambitious goal (i.e. implementing quantum information). 



1.1.1 The quantum bit. Classical information is based on binary logic, in 
which information is encoded in a series of bits (binary digits) that can assume 
only two values, or 1. A typical example is a switch, with its two possible 
states "on" and "off". All classical logical operations can be implemented as 
algorithms based on one- and two-bit operations, the so-called gates. 

The building block of quantum information is the quantum bit, or qubit. 
Using the Dirac notation, the two states that characterize the qubit are |0) 
and |1) and they represent the quantum counterpart of the classical and 1. 
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The most important property of a quantum bit is the possibiHty to be in a 
coherent superposition state 

|V) = «|0)+/3|1), (1) 

with a and /3 complex numbers, characterized by a relative phase and by 
|ap + |/3p = 1. According the the postulates of quantum mechanics, |ap 
represents the probability for the qubit to be in the state |0), whereas 
represents the probability to be in This means that if we prepare many 
copies of the same system in the state a measurement of the state of 
the qubit will produce the outcome with rate lap, and the outcome 1 with 
rate In which case the two logical states are the spin up | |) and the 

spin down | |). The two states |0) and |1) form a basis of the Hilbert space 
H. = span{|0), |1)} of the qubit. 

A good example of a qubit is the spin 1/2. In order to explain the necessity to 
use complex number a and /3 to characterize the state of the qubit, we describe 
an interference procedure for a spin 1/2 particle. Suppose we prepare the spin 
in the state |?/'o) = | T), that is with probability 1 to find it parallel with 
respect to a certain direction z in the space, that we choose as quantization 
axis. We then rotate the spin by an angle 7r/2 about an axis perpendicular to 
z, i.e. the y axis. The result is the state 

|^i)=e-^5'^1T) = -^(|T) + U)). (2) 

We then let the spin cross a region in which there is a magnetic field that 
points in the positive z direction, B = {0,0, B). Due to the presence of the 
magnetic field, the two states | j) and | j) accumulate a relative phase 2(f, 
that depends on the magnitude of the magnetic field and the time t spent in 
the region with the B field, and that for simplicity we leave unspecified. Up to 
an overall phase, the state of the system that comes out from the region with 
a magnetic field is given by 

1^2) = ^ (I T) + e^^n i)) . (3) 

Now, we again rotate the spin of 7r/2 about the y direction, and obtain 

IV3) =e-^?'^«|V2) =e^'^[cosM| T> + ismM| i)] . (4) 

If we now measure the state of the spin, we obtain | |) with probability cos^((/3) 
and I I) with probability sin'^ (ip). We clearly see, now, that the relative phase 
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can really affect the state of a quantum system. This procedure is known 
as Ramsey interference [1] and it is used in experiments to detect coherent 
oscillations in the transverse spin component. 

1.1.2 One qubit as environment. Decoherence is a consequence of the in- 
teraction of the qubit with the surrounding environment. As an instructive 
example we consider the case in which the environment is constituted by an- 
other qubit. For the Hamiltonian describing the interaction between the two 
qubits we choose {h = 1) 

H = ^al^al, (5) 

where the operator cif ®(t| is a two-qubit operator, given by the tensor product 
of two single-qubit operators, and it acts in the tensor product space TC = 
Til (8) 7^2- We note that our general argument does not depend on the specific 
form of H, as long as it describes an interaction between the qubits. As the 
initial state for the two qubit system we choose a product state |+)i |+)2, 
where the single qubit state is |+) = (|0) + |l))/\/2, written in the basis 
diagonal with respect to a^, az\0) = |0), and cr^ |1) = — We let the system 
evolve according to the unitary evolution generated by the Hamiltonian Eq. ([5]) 
for a time t, after which we perform a trace operation on the second qubit and 
have a look how the state of the first qubit evolved during the time t in which 
it has interacted with the second qubit. We re-write the initial state of the 
first qubit as a pure state density matrix, pi = |+)i(+|. In the {|0), |1)} basis 
it is found to be 

m(o)=i(;;). (6) 

The state of the two qubit system after a time t is given hy \ip{t)) = f/(t)| + +), 
with U(t) = exp(— i Jtcrf (t|/2). After some algebra the state of the first qubit 
at time t is given by 

= Tr,||^(t)>{V,(t)II = i ™(f /^)) . (7) 

The diagonal element of the first qubit density matrix is left unchanged by the 
interaction with the second qubit, whereas the off-diagonal elements change 
in time. The coherence of a state is encoded in the off-diagonal element of the 
density matrix. After a time t = n/J the coherence is completely lost (full 
decoherence). However, due to the smallness of the environment considered. 
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the first qubit periodically recovers its original state. It is therefore clear that 
the interaction with the environment strongly affects the qubit coherence. 

1.2 Quantum open systems 

According to the axioms of quantum mechanics the dynamics of a closed con- 
servative system is described as a unitary time evolution. In such a picture the 
system is considered to be decoupled from the surrounding environment, that 
does not influence at all the dynamics of the closed system. Strictly speak- 
ing this is never the case. However, under certain conditions the coupling to 
the environment can be considered to be weak, and to a good approximation 
neglected. 

In condensed phases, the coupling to the environment can be relatively 
strong, and the system under consideration cannot be separated from its sur- 
rounding. However, often a rather complex physical situation can be modelled 
by a system that consists of few dynamical variables in contact with a huge 
environment, constituted by a very large or even infinite number of degrees of 
freedom. In this case the small relevant system alone has to be described as 
an open system. 

Generally, an open system is a quantum system S which is coupled to an 
other quantum system B called the environment. It can therefore be thought 
to be a subsystem of the combined system S + B, which in turn in many 
cases is considered to be a closed system, governed by Hamiltonian dynamics. 
The system S will in turn change according to its internal dynamics, and 
as a consequence of the interaction with the environment. Certain system- 
environment correlations will be established between the two parts and, as a 
consequence, the dynamics of a quantum open system cannot, in general, be 
described in terms of a unitary time evolution. 

Denote Hs the Hilbert space of the system S, and Hb the Hilbert space 
of system B. The dynamics of the combined system S + B takes place in the 
Hilbert space given by the tensor product space Hsb = T~(-s 'Si'Hb- The total 
Hamiltonian can be chosen to have the general form 

H = Hs + Hb + Hi, (8) 

where Hs describes the evolution of the system S alone, Hb is the free Hamil- 
tonian of the environment B, and Hj describes the interaction between the 
system and the environment. Usually when speaking about the environment 
of the system S, the term reservoir may appear, that refers to an environment 
with an infinite number of degrees of freedom, such that the frequency modes 
associated with it form a continuum spectrum. Occasionally the term heat bath 
or simply bath refers to a reservoir which is in thermal equilibrium. 
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The presence of an environment is meant to model the communication of 
the open system with the external world. However the attention is focused 
on the subsystem under study S, and all observations of interest refer to the 
subsystem S. Formally this means that all observables of interest act on the 
Hilbert space Hs- Denoting the state of the total system by p, the expectation 
values of all observables may be written as 



{0)=Trs[Ops], ps = Tvb[p], 



(9) 



where O is the Hermitian operator describing the observable, ps is the reduced 
density matrix of the open system S, and Trg^g) denotes a partial trace on 
the system S{B). 

All informations that describe the open system S are contained in the re- 
duced density matrix ps- Since the total system evolves unitarily in time, ps{t) 
is obtained as partial trace over the environment B of p{t), 



ps{t) = TTB Uit,to)psito)U\t,to) 



(10) 



where U{t, to) is the unitary evolution operator of the total system. The equa- 
tion of motion for the open system reduced density matrix ps{t) is 



dt 



ps{t) = -iTxB [H{t),ps{t)\. 



(11) 



1.3 Generalized master equation 

In many cases it is useful to model the dynamics of an open system by means of 
an appropriate equation of motion for its density matrix, the so called quantum 
master equation. The evolution in time of the total system p is governed by 
the well known Liouville equation of motion 

p = -i[H{t),p{t)] = -iCp{t), (12) 

where the second equality defines the Liouville operator C As the Hamiltonian 
can be divided into three terms that describe the dynamics of the two systems 
alone, Hg and Hb, and a interaction between the two parts, i?/, the Liouville 
operator can be written as the sum of three contributions 

C = Cs + Cb + Ci. (13) 

Without going into details that arc beyond the scope of this review, we just 
mention the fact that the Liouvillian is a superoperator, that maps operators 
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into operators. The initial state for the combined system S + B can typically 
be chosen to be a product state, p{0) = ps{0) Pb- 

We have already introduced the reduced density matrix ps of the open 
subsystem S. It can be formally obtained from the density matrix of the total 
system p by means of a projection operation, that contains a partial trace over 
the system B, 

PS = Vp = Ttb[p\® PB- (14) 

Here, pB is a fixed density matrix for the environment. Mapping operators into 
operators, the projector V is also a superoperator. We may thus decompose 
the p as 

p{t) = ps{t) + {l-V)p{t), V^=V. (15) 

Substituting this decomposition in the Liouville equation of motion for the 
total system Eq. (jl2p . choosing the projector in such a way that the inhomo- 
geneous term that depends on the initial state can be disregarded, and using 
that the operator V defined in Eq. (|14p commutes with the Liouvillian of the 
open system Cs-, after some algebra, the equation of motion for the reduced 
density matrix ps{t) can be cast in the form of an exact generalized master 
equation, the Nakajima-Zwanzig equation [2,3] 

ps{t) = -iCsPsit) + [ dt'j:{t - t')ps{t'), (16) 
Jo 

S(t)/55 = -iTTB [Cie^^-'^^'^'Cips ® pb] , (17) 

where S(t) is the self-energy superoperator. The first term describes the re- 
versible evolution of open system S, while the second term produces irre- 
versibility. 



1.3.1 Born approximation. The generalized master equation Eq. (|16p is 
a formally exact and closed description of the dynamics of the state of the 
system ps, but it is very complicated from a mathematical point of view and 
rather unpractical. Usually, in order to handle it and some approximation 
are made. In fact, the kernel of Eq. (I16p contains all powers of £/, and the 
dynamics of ps at time t depends on the whole history of the density matrix. 
If the coupling between system and reservoir is weak, i.e. ^ \\Cs + ^_b||, 
the exponential can be expanded in power of £/ in a perturbative way. In 
lowest order Born approximation, the interaction Liouvillian is disregarded in 
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the exponent and Ci is retained only to second order 

t{t)ps = -iTiB [£/e(i-^)(^-+^-)*£/p5 ® Pb] ■ (18) 

The apphcabihty of the master equation in the Born approximation is strictly 
restricted to those cases in which the coupling between system and environ- 
ment is weak, with decoherence and relaxation times large compared to the 
relevant time scales of the reversible dynamics. 

1.4 Quantum Markov process 

The master equation in the Born approximation Eq. (jlSp . though it is much 
simpler than the exact Nakajima-Zwanzig equation (|16p. is still an integro- 
differential equation that is very difficult to handle. Assuming that the tem- 
poral correlations in the bath are short lived and typically lead to exponential 
decay of the coherence and populations, the master equation in the Born 
approximation Eq. (jlSp can be further simplified. In the Born-Markov ap- 
proximation the master equation for the reduced density matrix of system S 
assumes the form 

ps{t) = -iCspsit) + t'^{t)ps{t), (19) 

t^{t) = -i f' dt't{tY''^'. (20) 
Jo 

In an eigenstate basis of Hs, the master equation in the Born-Markov ap- 
proximation can be written as the so called Redfield equation [4-6] 

Pnm = —i^nmPnm{t) — ^ RnmklPkl{t), (21) 

k,l 

where pnm = {fAps\m)^ ujnm = ^n — ^mi and we have introduced the Redfield 
tensor 

Rnn^ki = / dtTTB [{n\ [Hf\t), [Hf\0), \k{t)){l{t)\pB]] \m)] , (22) 
Jo 

where we have used the interaction picture Hamiltonian and the system eigen- 
states in the interaction picture 

Hf\t) = e^(-f^^+-f^«)*i/^e-^(^^+^«)*, \k{t)) = e'^'^\k) = e^'^'=*|A;). (23) 



September 26, 2008 21:16 Advances in Physics DecoSSQ-rcview 



The first term of Eq. (j2ip represents the reversible motion in terms of the 
transition frequencies LOnm, while the second term describes relaxation. The 
Redfield tensor can be expressed as 

Rnmkl = ^nm ^ ^nrlk + ^n-k ^ ^Wrm ~ ^\mnk ~ ^^Imnk^ (^^) 
r r 

in terms of rates given by the golden rule expression 

^l^nk = / dte-^^-^TrB Hl{t)lraHmnkPB , (25) 
Jo L J 

TztL = r dte-'^--'TTB \Hj{0)lraHl{t)nkPB] , (26) 
Jo L ^ 

with Hi{t)im = (n|e^*^-i?7e-^*^-|m), and {ri^ik)* = rJ-J^^. 

We have already pointed out that the dynamics of an open system cannot 
be described as a unitary evolution. However, the mapping describing the 
evolution is required to be completely positive [7], implying p — > Yln^npOt, 
where {On} is a set of linear operators on the reduced state space that satisfy 
OnOn = 1, such to preserve the trace of p. In the framework of Lindblad 
theory [7], the master equation can be cast in the form 

ms = -i[Hs,ps{t)] + \ ^ {[L,p5(t),4] + [Lj.Psm]]} . (27) 

i 

The Lindblad operators Lj describe the effect of the environment in the Born- 
Markov approximation. 



1.4.1 Two level systems and Block equations. The aim of this review is to 
provide a overview on the mechanisms that affect qubit dynamics and induce 
decoherence in solid state realizations of qubits. Therefore we concentrate on 
two level systems and their coupling to the surrounding environment. 

The density operator of a two state system is a two dimensional positive 
Hermitian operator with trace one. It can thus be expressed in terms of a basis 
of Hermitian operators given by the three Pauli operators cr = {ax,cry, az) and 
the 2x2 identity, 



P= ^(1+P-cr), 



/ Poi + Pio \ 
p = Tr[pcr] = i{poi - pio) I 

V Poo - Pu J 



(28) 
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The vector p is known as the Bloch vector, and for a spin- 1/2 object it repre- 
sents the expectation values of the spin components p/2 = (S) = Tr[Sp], where 
S = cr/2, with cjz diagonal in the |0) |1) basis, cr^lO) = |0) and cr^jl) = — 11). 
Combining the last equation with the Redfield equation (j2ip in the case that 
n,m,k,l = 0, 1, the master equation within the Born-Markov approximation 
for the density matrix of a two level system can be expressed as a first or- 
der time differential equation for the expectation value of the spin component 

{S) = i{S,),{Sy),{Sz)), 

(S)=^x (S)-ii(S) + (So), (29) 

with LU = (0, 0,u;oi). In case of a spin 1/2 particle in a magnetic field defin- 
ing the z direction, wqi represents the Zeeman splitting. The inhomogeneous 
term (Sq) and the relaxation matrix R depend on the rates Eqs. (I25p . (I26p . If 
(Joi S> Rnmki, it is possible to make a secular approximation, retaining only 
terms Rnmki with n — m = A; — /, [8], such that the Redfield tensor can be 
approximated by the diagonal form 

/T^-i \ 
i? Tg"^ , (30) 

V TfV 

where the relaxation time Ti and the decoherence time T2 are given by 



^ = 2Re(rSo + rS;^i), (31) 
111 

— = + — , 32 



with the pure dephasing time T,^, given by 

— = Re(rQQQQ + r[-^^^-^ — 2rQQ^^-^). (33) 



For a system-environment coupling given by a simple bilinear form Hi = 
Os ® Xb, with Os an operator acting in the system space TLs, and Xb an 
operator acting in the environment space "Hb, the relaxation and dephasing 
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times Ti and can be written as 



^ = 4KO|05|l)|V(a;oi)coth^, (34) 



2kBT, (35) 



where the spectral density J{uj) is the Fourier transform of the environment 
time correlator 

/oo 
dtTrB[XBXB{t)pB]e-'^'. (36) 
-oo 



The first term in Eq. (j29p produces a rotation of the Bloch vector along the z 
direction. If i? = we have the classical picture of a magnetic moment process- 
ing along the externally applied magnetic field. The second term proportional 
to R describes an exponential damping of the component of the Bloch vec- 
tor. Ti describes the decay of the longitudinal component of the Bloch vector, 
while T2 describes the decay of the transverse component. 

We remark that the Markovian results Eq. (j32p satisfy the expected funda- 
mental Korringa relation [9]. 



1.5 Spin-boson model 

Here we describe a simple model to treat the dynamics of a two-level system 
in contact with a reservoir. We consider a generic two-level system described 
by the Hamiltonian 

■Hs = ^a, + '-a,. (37) 

In order to include the effect of dissipation in the quantum formalism, it is cus- 
tomary to follow the Caldeira-Leggett [10-12] approach. A bath of harmonic 
oscillators at thermal equilibrium at temperature T is introduced to describe 
the degrees of freedom of the environment. The system-l-bath Hamiltonian is 

H = Hs + Hb + Hsb, (38) 
HB = \Y.u:^{bib^ + ^, (39) 

HsB = O5 ® = ^ c« (he, + , (40) 
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where TCs is the quantized Hamiltonian of the system Eq. (j37p . TCb is the 
bath Hamiltonian, described by independent bosonic degrees of freedom with 
frequencies uja- The couphng between the system and the bath degrees of 
freedom is described by Hsb, where Os = cFz, Xb = Yla + ^q)) Cq, 
are couphng parameters. 

A rigorous treatment of the spin-boson model in the Born approximation 
without making use of the Markov approximation is presented in [13, 14]. The 
eigenstates of the Hamiltonian Eq. (j37p are 

|0) = ^ (./l + ^|+) + ./l-— I-)) , (41) 
V2 VV ^01 V "^01 / 



1 

72 V 



VV ^01 V "^01 / 



(42) 



where |ib) are eigenstates oi Uz, o"^|ib) = ±|ib), and ojqi = VA^ + e^. The 
initial state of a system of quantum harmonic oscillators in thermal equilibrium 
is 

PB = exp(-/??^B), Zb = Trexp(-/3HB). (43) 

All the informations of the bath, such as the bath frequencies and the 
coupling parameters Cq appearing in the Hamiltonian, are contained in the 
spectral density J{oj) of the system-bath coupling, 

a 

Here we limit our attention to the Markovian case, and make use of the general 
Redfield theory described in the previous section. From the formula Eqs. (j34p . 
()35p . the relaxation and dephasing rates take the form 

2 

^(c^oOcoth^, (45) 





i- 


Ti 


V'^Ol 








V^^Ol 



' J(^) 



UJ 



2kBT. (46) 



1.6 Spin qubits 

Per antonomasia, the two-state system that nature provides us with is the in- 
trinsic angular momentum of the electron: the spin 1/2. It is therefore natural 
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to choose the electron spin as the two-state system that encodes the qubit. 
The spin of the electron can have much longer decoherence time than the 
charge degrees of freedom. Nevertheless, isolating the spin degree of freedom 
of an electron to a degree required for quantum computation is not at all an 
easy task. Moreover, in order to be used for quantum computational purposes, 
electron spin-based qubits must be designed as scalable devices that can be 
externally controlled, coupled, manipulated, and read-out, i.e. they must sat- 
isfy the DiVincenzo criteria [15]. A successful and promising device for the 
physical implementation of electron spin-based qubits is the semiconductor 
quantum dot [16]. 

1.7 Semiconductor quantum dots 

The quantum dots owe their name to the zero-dimensional character of such 
devices. Can be considered as a quantum box that can be filled with electrons 
(or holes) which occupy the available discretized states of the system. The 
electrons can tunnel on and off the dot, which is coupled to large reservoir 
via tunnel barriers. The height of the barriers, and consequently the rates for 
tunneling through the barriers on and off the dot, can be controlled via the 
application of gate electrodes. Electrostatic gates can also be used to tune the 
electrostatic potential of the dot with respect to the reservoirs, such that the 
ladder of energy levels in the dot can be shifted up or down with respect to the 
energy of the reservoir. External bias voltages can be applied and transport 
properties can be measured. 

Quantum dots are basically characterized by the quantized level structure, 
for which they are considered as artificial atoms, and by the transport state 
of the dot, that can be active or blocked, and depends on the combination of 
bias and gate voltages applied. In fact, the Coulomb repulsion between the 
electrons in the dot determines an energy cost for adding an extra electron in 
the dot. At low temperatures, the tunneling of electrons on and off the dot can 
be drastically suppressed, and the dot is in the so called Coulomb blockade. 

Many kinds of quantum dots have been realized so far. Here we focus the 
attention on lateral III-V semiconductor quantum dots, as those in Fig. [TJ 
These devices are fabricated from heterostructures of GaAs and AlGaAs grown 
by molecular beam epitaxy. The energy potential along the growth direction of 
such a structure has a minimum at the interface of the two layers, which is also 
asymmetric with respect to the growth direction. Free electrons are introduced 
by doping the AlGaAs layer with Si, which accumulate at the GaAs/ AlGaAs 
interface, deep down in the minimum of the vertical potential, that provides 
strong confinement of the electrons along the growth direction. At the same 
time, the electrons are free to move along the interface, where they form a 
two dimensional electron gas (2DEG), that can have a high mobility and a 
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Figure 1. Scanning electron micrograph of a double quantum dot. [17] (With permission from the 

Authors). 



relatively low electron density (typically 10^ — 10^ cm^/Vs and ~ 10^^ m~^). 
The low density results in a relatively long Fermi wavelength 40 nm) and a 
large screening length, such that via application of an electric field, obtained 
through metal gate electrodes on top of the heterostructure negatively charged, 
the 2DEG can be locally depleted. Therefore, by suitable designing the gate 
structure it is possible to isolate small islands of the 2DEG, thus creating a 
dot. When the lateral size of the dot is compared to the Fermi wavelength, 
the energy level structure of the dot becomes discretized, and at temperatures 
down to tens of mK, the energy separation of the levels becomes much higher 
than the temperature, such that quantum phenomena start to play a significant 
role. 



1.8 Spin relaxation and spin dephasing mechanisms in quantum dots 

The electron spin in semiconductor quantum dots can be isolated and con- 
trolled with a high accuracy, but it still suffers from decoherence due to the 
unavoidable coupling with the surrounding environment. In order to imple- 
ment quantum computation algorithms with electron spin-based qubits in 
semiconductor quantum dots, it is necessary to engineer the devices in such a 
way to preserve the coherence of the electron spin states for sufficiently long 
time scales. Besides the fundamental interest, it is therefore important to the- 
oretically understand which sources of dissipation and decoherence affect the 
electron spin in quantum dots, and to find ways to reduce their infiuence on 
the spin-qubit dynamics as much as possible. 

Two kinds of environment turn out to mainly afi"cct the dynamics of an 
electron spin in a quantum dot, the phonons in the lattice, and the spins of 
atomic nuclei in the quantum dot. 
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Phonon-induced relaxation in semiconductor quantum dots has attracted 
some attention from a theoretical point of view for the first time in [18, 19]. 
The lattice phonons do not couple directly to the spin degree of freedom. 
However, even without the application of external electric fields, the breaking 
of inversion symmetry in GaAs gives rise to spin-orbit interaction, that cou- 
ples the spin and the orbital degrees of freedom. The latter, being coupled to 
the phonons, provide an indirect coupling between the electron spin and the 
phonons, that constitute a large dissipative bosonic reservoir and provide a 
source of decoherence and relaxation. Short time correlations in the phonon 
bath induce a Markovian dynamics of the electron spin, with well defined re- 
laxation and decoherence time Ti and T2. It turns out that effectively the 
phonon-induced pure dephasing time of an electron spin in a quantum 
dot in the presence of a magnetic field diverges. In the Bloch picture, pure 
dephasing arise from longitudinal fluctuations of the magnetic field, while a 
perturbative treatment of the spin-orbit interaction gives rise, within first or- 
der, to a fluctuating magnetic field perpendicular to the applied one. As a 
consequence the decoherence time T2 is limited only by its upper bound Ti, 
T2 = 2Ti. In turn the relaxation time Ti shows a strong dependence on the 
magnetic field. Ti oc B^, that has been confirmed experimentally, where a very 
long relaxation time up to Ti « 1 s has been measured for a magnetic field of 
B = 1T [21]. 

Hyperfine interaction was for the first time taken into consideration as a 
source of decoherence for an electron spin confined in a quantum dot in [22]. 
In GaAs there are 5 x 10^^ atoms in 1 cm^. Therefore, the linear extension 
of a typical GaAs quantum dot, that is of order of the Fermi wavelength 
~ 40 nm, encompasses roughly 200 atoms, from which it can be estimated 
that the wavefunction of an electron in a GaAs quantum dot overlaps with 

10^ nuclei. The electron spin and the nuclear spins in the dot couple via the 
Fermi contact hyperfine interaction, that creates entanglement between them 
and strongly affect the electron spin dynamics. It turns out that long time 
correlations in the nuclear spin system induce a non-Markovian dynamics of 
the electron spin, with non exponential decay in time of the expectation values 
of the electron spin components. In a large applied magnetic field B, the 
dynamics in the nuclear field due to the hyperfine interaction can be treated 
perturbatively and it turns out that flip-flop dynamics starts to affect the 
nuclear field in a time that scales like the number of nuclear spins, oc A^. For 
shorter times the nuclear field is static and the transverse component shows 
a Gaussian decay, that is due to the statistical distribution of nuclear spin 
states. 

We remark that the phonon induced relaxation rate of the electron spin is 
enhanced by an applied magnetic field, whereas the influence of the hyperfine 
interaction is reduced by a large Zeeman splitting. 
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1.9 Hyperfine-induced decoherence in spin qubits 

1.10 Hyperfine interaction 

The spin of an electron and the atomic nuclear spin can interact through the 
hyperfine Fermi contact interaction, a spin-spin interaction that takes place 
when the electron and the nucleus occupy the same position in space, from 
which the term "contact". The origin of the hyperfine coupling can be un- 
derstood considering the electromagnetic interaction of an electron with the 
magnetic field produced by a nucleus. Without loss of generality the mag- 
netic properties of a nucleus can be described as those of a magnetic dipole 
/^Af = (J-nKI, where fiN is the nuclear magneton, and I is the nuclear spin op- 
erator. The interaction of a nuclear dipole with the electronic shell gives 
a rather small effect, and can be treated using a perturbative method. In the 
non-relativistic Pauli description of the electron, the Hamiltonian of an elec- 
tron in a magnetic field B = V x A produced by a vector potential A is given 
by 

W = ^(p + ^a)' + 2^bS-(Vx A), (47) 

where S is the electron spin operator. The vector potential produced by a 
magnetic dipole fj, at position r is, according to classical electromagnetism, 
A = (/I X r)/r^ = V X (/i/r). Neglecting the term quadratic in the vector 
potential and replacing KL = r xp for the electron orbital momentum operator, 
the Hamiltonian Eq. (j47p can be written as 

n = 2^B^^^ + 2/iB(S • V)(/i • V)- - 2^b(S • /x)V2-. (48) 

The magnetic interaction of the nuclear spin and the electron spin is contained 
in the second and the third term of Eq. (|48p , and it is obtained after integration 
over the orbital degrees of freedom, i.e. it has to be understood as applied to an 
electron orbital state ^/'ei(r). For r 7^ 0, the terms involving the electron spin S 
in Eq. (j48p behave regularly, the last term vanishes identically, while the second 
term produces a usual dipole-dipole interaction 2/i5[3(S-r)(/x-r)/r^ — S-///r^]. 
The case r = needs to be treated more carefully. It can be shown [9] that 
the dominant contribution of the spin dependent part of Eq. (j48p reduces to 
(167r/3)^B(S • /i)(5(r), and once applied on the electron orbital wave function 
is given by 



^hy = Yvr;^B|V'el(0)|^S • /X, 



(49) 
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Figure 2. Schematic representation of the orbital wave function of an electron in a quantum dot. 
Due to the spatial extension of the wave function, the electron spin (big blue arrow) couples to 
many nuclear spins (small red arrows). 



which is finite for s electrons and zero for others. The Hamiltonian for the 
magnetic interaction of the electron with the nucleus can be written as 



(50) 



1.10.1 Hyperfine interaction in semiconductor quantum dots. In a semi- 
conductor quantum dot an electron is confined in a 2D region of space whose 
linear extension is of the order of the Fermi wavelength, for GaAs about 
~ lOOnm, that is much larger than the typical lattice spacing of the crys- 
tal (~A). As a result a discretization of the energy levels in the dot appears, 
with an orbital level spacing that, for lateral quantum dot containing single 
electrons, is much greater than the typical energy scale of the hyperfine interac- 
tion. As opposed to the case of single atoms, the electron orbital wave function 
in a quantum dot extends over a region much larger than the lattice size, such 
that the electron spin couples to many nuclear spins, as schematically shown 
in Fig. [2j The effective hyperfine Hamiltonian describing the interaction of a 
single electron with the nuclei in the dot can be written as 

N N 

n = S-h, h = AvoY,\M^k)\% = Y.^klk, (51) 

fc=l k=l 

where S is the electron spin operator, h is the so-called Overhauser field, 
given by the sum of all the Ik nuclear spin operators, weighted by the position 
dependent coupling strength = wo^|^o(rfc)|, where the square modulus 
of the electron envelope wave function at the kth. lattice site. Typically the 
electron can be assumed to be in the quantum dot orbital ground state, vq 
is the volume of the crystal unit cell containing one nuclear spin and A = 
IdirfiBfJ'N^/'i is the contact hyperfine coupling strength. In GaAs the nuclear 
spin is / = 3/2 and an estimate of the interaction strength, weighted by the 
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abundances of the three isotopes naturally present (^^Ga, ^^Ga, and '''^Ga), 
yields A ^ 90|UeV. 

The inhomogeneity of the electron wave function results in a nonuniform 
hyperfine coupling strength Ak, that depends on the probability to find the 
electron in the nuclear lattice site k, resulting in a subtle and complex many- 
body quantum mechanical behavior, with the nuclear spin affecting the elec- 
tron spin time evolution, and the electron spin acting back on the dynamics 
of each of the nuclei in turn. 

From the point of view of the electron spin, entanglement with the degrees 
of freedom of nuclear spins arising from the hyperfine coupling constitutes a 
decoherence mechanism. 

Since the Overhauser field h appearing in Eq. (fSTIl is composed by the sum of 
a large number of spins, it is natural to question whether the nuclear field can 
be approximated to a classical object and to which extent this approximation 
gives correct results. In a relatively recent work by Coish et al. [23] it has 
been theoretically shown that, for the special case of a uniform hyperfine 
coupling constants A}^ = A/N, arising from a constant wave function in the 
dot ipo = 1/\/vqN, the dynamics obtained in the mean field approximation and 
the quantum evolution show agreement up to the transverse-spin correlation 
time Tc, which diverges in the zero external magnetic field case (unphysical 
result due to the assumption of constant coupling), but that saturates to a 
finite value in case of a finite external magnetic field. 



1.10.2 Fluctuation timescales of the nuclear field. The nuclear field is 
quantum many-body interacting spin system whose field orientation and mag- 
nitude change over time. This change is due to the combined effect of the 
inter-nuclear dipole-dipole interaction and the hyperfine interaction between 
electron and the nuclei. The dipolar interaction does not conserve the total 
nuclear spin and thus can be responsible for changes in the nuclear spin config- 
uration. Those changes, combined with the spatial variation of the hyperfine 
coupling constant, lead to a different value of the nuclear field seen by the elec- 
tron spin and thus to its decoherence. Here we briefly outline the timescales 
in which those mechanisms take place, in order of decreasing timescales. 

The strength of the effective magnetic dipole-dipole interaction between 
neighboring nuclei in GaAs is directly given by the width of the nuclear mag- 
netic resonance (NMR) line to be ~ (lOO/is)"^ [24] and its inverse can be taken 
as an estimate for the timescale in which a change in the nuclear configuration 
due to dipolar interaction takes place, i.e. Td-d ~ which is just the pe- 

riod of precession of a nuclear spin in the local magnetic field generated by its 
neighbors. This timescale is so long that a great number of other decoherence 
mechanisms start to play a significant role before nuclear dipole effects start 
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to matter. 

Besides spin diffusion driven by nuclear dipole-dipole interaction, the nu- 
clear field can undergo a change due to the flip-flop term in the Hamiltonian 
Eq. (jSip . In a large external field B, the flip-flop term can be treated within 
the framework of perturbation theory, as it will be explained in the next sec- 
tion. We anticipate here that the hyperfine mediated dynamics in the nuclear 
field has a timescale given by cx A/N. This means that up to this timescale 
the nuclear field can be considered as static. 



1.11 Decoherence due to hyperfine-induced electron spin dynamics 

An early treatment of the hyperfine interaction as a decoherence mechanism 
for single electron spins confined in quantum dots whose carried out in [22]. 
There, a second order time-dependent perturbation expansion of the hyperfine 
interaction in a magnetic field was performed with respect to the flip-flop 
transverse term A{hj^S- + /i_S'+)/2 for a constant hyperfine coupling A and 
a long-time longitudinal spin-flip probability ~ is obtained, where p 

is the nuclear spin polarization. As a result beside a large external magnetic 
field, a large polarization p and a large number of nuclei in the dot would 
suppress the spin-flip probability. 

The first signature of the non-Markovian behavior of the nuclear spin bath 
appeared in [26], where an exact solution for the fully polarized case p = 1 
was provided. In that case the decoherence is due to a non-uniform hyperfine 
coupling that depends on the probability for the electron to be located at dif- 
ferent nuclear sites. A remarkable feature of the non-Markovian behavior is the 
long-time power law decay of the electron spin correlator, ~ in strong 

Zeeman field, according to which the longitudinal electron spin component 
decays of a fraction of ~ in a time ~ N/A. 



1.11.1 Single- electron spin decoherence in large Zeeman splitting. A de- 
tailed and comprehensive treatment of the hyperfine interaction [28] provides 
an analytical result for the electron spin dynamics for arbitrary nuclear spin 
/ and nuclear polarization p. A generalized master equation (GME) approach 
allows a treatment of the transverse electron spin-nuclear spin flip-flop terms 
in the Hamiltonian with an external field in a well controlled perturbative way. 
An expansion of the self-energy in the exact Nakajima-Zwanzig GME shows 
a rich electron spin dynamics, with exponential and non exponential decaying 
contributions and undamped oscillations. The form of the decay of the trans- 
verse and longitudinal electron spin component is obtained in high magnetic 
field up to forth order in perturbation theory. 
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The hyperfine Hamiltonian in an external magnetic field is 

H = hS, + en J, + h • S, (52) 

where b = g* hbBz {enz = 9WnBz) is the electron (nuclear) Zeeman splitting 
in a magnetic field defining the z-axis i?2, g* (gj) the effective electron (nuclear) 
^-factor, and hb (fJ-N) the Bohr (nuclear) magneton. In the rotating frame with 
respect to the nuclear Zeeman term the Hamiltonian can be separated into a 
longitudinal (unperturbed) and transverse (perturbation) term, 




In absence of V, {Sz)t is constant, since [Ho,Sz] = 0, but the transverse 
component {S±)t evolves in time in a non trivial way. For a large number of 
nuclear spins ~ 10^ (GaAs dot) a direct application of the central limit 
theorem gives a Gaussian distribution for the eigenvalues of hz with mean 
ho = {hz) and variance a ~ A/\/l^. The transverse correlator for an initial 
state given by the product state of the initial electron spin state psiO) and 
incoherent Gaussian distributed nuclear mixture state is 



{S+)t « (5+)oexp[-tV2r2 + ^(6 + /io)t], r = - = ^Jj^- (54) 

Choosing as nuclear initial state the pure state 1-0/(0)) = J|^-(^l + p\ |j) + 
g*03^1 — p\ [j))/\/2, for a certain polarization p, the same result Eq. (j54p with 
ho = pN comes out. 

The reason for this decay lies in the choice of the initial nuclear state con- 
taining many hz eigenstates and can also be obtained choosing the nuclear 
field in a hz eigenstate, but with the electron spin in a transverse initial state. 
This decay is reversible and can be removed with a standard spin echo tech- 
nique [29,30]. Such an experiment therefore reveals only the decay due to the 
transverse flip-flop term V Eq. ()52p . A procedure more suitable for a quan- 
tum computation algorithm would be a strong Von Neumann measurement of 
the nuclear field that would then prepare a hz eigenstate, leading to simple 
precession with no decay [31,32]. 

Analysis of the GME in the Born approximation for a very high magnetic 
field (6 S> N) provides an asymptotic form to leading orders in ~ 1/uJn = 
yib + K), 

{S+)t « (t) + crf'it), {Sz)t ^ {Sz)oo + af^it), (55) 
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Figure 3. Decay of the driven Rabi oscillation in ESR showing a power law decay oc l/\/t and a 
universal phase shift of 7r/4, as given by Eq. JSTj, valid for t 2> max(l/(T, l/feac,t'ac/2(T ). For the 
plot the value b = a = 0.4 Ghz have been chosen. 

where {S,)oc oc {S,)o, oc {S+)oe'^-\ and ajf^(t) oc for a 

parabolic confinement in the dot, with (5 <C 1. Even for a hz eigenstate, for 
which no decay is expected in zeroth order in the transverse electron nuclear- 
spin flip-flop interaction, a long time irreversible decay takes place, that is due 
to the spatial variation of the hyperfine coupling constant. 

1.11.2 Single-spin ESR: universal phase shift and power law decay. Here 
we describe the situation in which the electron spin is coherently driven via 
pulsed magnetic resonance, while coupled to a nuclear long-time correlated 
spin bath. Recent remarkable experimental results [30] show a coherent elec- 
tron spin oscillation, even for a Rabi period much longer than T| = 10 — 20 ns. 
A non-exponential decay of the Rabi oscillations is observed, obeying a power 
law decay with the appearance of a universal phase shift. 

Consider a quantum dot in a time independent magnetic field defining the 
z direction. In addition, an oscillating magnetic field is applied in the plane, 
along the x direction. For a large number of nuclei in the dot {N ~ 10^ in GaAs 
dots) the field is Gaussian distributed, with mean /iq and variance a [25,26, 
28]. In the case of strong external field {b ^ a, with b = g* /j^bBz), neglecting 
transverse electron-nuclear spin fiip-flop terms in the hyperfine interaction, the 
Hamiltonian is (h = 1) 

TC = {b + hz)Sz + ba,cCos{ujt)S::,, (56) 

where 6ac = 5*/^B-Bac, uj and i?ac being the frequency and amplitude of the 
ESR driving field, is considered static (justified for t < Ifjs), and the 
assumption u = 6 + /iq is made. In the rotating wave approximation (valid for 
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ip'Ac/b)'^ 1), after averaging over the Gaussian distribution of hz, the decay 
of the driven Rabi oscihation is given by [30] 

for t » max(l/cr, l/6ac) &ac/2cr^), with C = 

exp(6^^/8(T^)erfc(6ac/\/8o-)\/27rfeac/8iT. The remarkable features appear- 
ing in the experiment arc the ~ l/v^ power law decay and the universal 
7r/4 phase shift. The reason for the appearance of these features is that the 
nuclear field hz does not change over a timescale much longer than the Rabi 
period. Since different values of hz determine different oscillation frequencies, 
an average over the distribution in give rise to a decay in the coherence 
of the driven electron spin, and the off-resonant contributions also determine 
the phase shift. The fact that coherent Rabi oscillations are visible even when 
the Rabi period is much longer than the transverse spin decay time r ~ 15ns 
has its origin in the fact that the power law decay sets in already after a short 
time l/c" ~ 15ns. 

In order to measure the electron spin state in the experiment [30], a spin- 
charge conversion technique is implemented by operating a double quantum 
dot in the spin blockade regime [33,34], in which the transport through the 
dots can occur only via transitions from spin states with one electron per 
dot, 1(1, 1)), to the singlet state in the right dot, |(0, 2)). The Pauli exclusion 
principle, that does not allows two electron with same spin state to occupy the 
same orbital, allows transport only for antiparallel spins. Transport of states 
spin-triplet states is therefore blocked. The oscillating transverse magnetic 
field rotates the spins, therefore unblocking an initial state with even parity 
spin state [34]. 

1.11.3 Single-triplet decoherence in a double quantum dot. An alternative 
way to implement a qubit with electron spin in quantum dots is to consider 
a double quantum dot with two spins, one per dot, and encode the qubit in 
the subspace with zero z-projection of the total spin S^^^ — Sf -\- S2 — 0. 
Advantages of this scheme is the possibility of reducing the hyperfine coupling 
in case of symmetric dots. At the same time additional decoherence due to the 
coupling to the orbital degree of freedom and leakage errors may appear. 

The cS'ective Hamiltonian for the one-electron-per-dot configuration can be 
written as 



Hdd = ezSz + h-S + Sh-6S + -S-S- J, 



(58) 
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where S = Si + S2, 5S = Si — S2, h = hi + h2, and 5h = hi — h2. Here J is the 
Heisenberg exchange couphng between the two electron spins. For definiteness 
we work in a regime of large Zeeman splitting due to an external magnetic 
field, £2 = g*^iB > max{((5h)i.ms, (h)rms}, where (C')rms = {i^i\0\ipi)^/'^ de- 
notes the root-mean-square expectation value of the operator O on the nu- 
clear state IV'/)- Requiring » J, where J is taken to be positive without 
loss of generality, the large Zeeman splitting condition renders the relevant 
spin Hamiltonian block diagonal with respect to the eigensubspaces of Sz- In 
the ^2 = subspace the spin Hamiltonian for the singlet \S) and 5^ = 
triplet I To), to zeroth order in the inverse Zeeman splitting l/e^, is given by 
T~(-o = (<^/2)S • S -|- Sh^5S^. The effective qubit Hamiltonian in terms of the 
vector consisting of Pauli matrices r = (r^,r^,r^), with the computational 
states |5) = —1) and |To) = 1) , has the form 

7^o = + + (59) 

A systematic treatment of the dynamics induced by the Hamiltonian Eq. (j59p 
can be found in [31,35]. The eigenstates of Ti.Q are given by a product state 
between a nuclear eigenstate |n) of dh^ and a superposition of \S) and \Tq), 
therefore Tio does not lead to any dynamics in the nuclear field. The correlator 
CtoS is defined as the probability to find the electron spins in the state [To) 
at time i > 0, provided the initial state {t = 0) was |V'(0)) = \S) \'>Pi), with 
IV'/) a superposition of 6h^ eigenstates, 

Cnsit) = 5^p/(n)Kn| ® {To\e-^'^°'\S) \n)f, (60) 

n 

where pi{n) diagonal matrix element of |'i/'/)(V'/| iii ths {jra)} basis. For a 
Gaussian distributed field dh^, with mean xq and variance do, the asymptotics 
of CtoS saturates to finite value that deviates from the semiclassical results 
[C™1(qq) ^ ;L/2] for J < xo [35] 

CT„5(oo)~H-Ki)'' (61) 
[ 2{J)\ ao<Cxo«J. 

At short times CTos{t) experiences a Gaussian decay on a timescale 
y\P + 4x^/Axoao, while in the case of strong coupling J S> max{Xo,cJo} 
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at long times t ^ J/^ctq a power law decay appears, [35] 

Cns{t) ~ Cto5(oo) - ^— ^ COS ^Jt + — J . (62) 

Those results show that the singlet-triplet correlator decays due to the quan- 
tum distribution of the nuclear spin system, even for a static system. For 
non zero exchange interaction J ^ the asymptotic behavior of the correlator 
CTosit) changes from a short time Gaussian behavior to a long time power-law 
(~ decay and acquires a universal phase shift which is 37r/2, consis- 

tent with experimental findings for the correlator Css{t) [36]. Qualitatively 
similar results appear when looking at the transverse correlator in the Sz = 
subspace, though one finds different decay power and different value of the 
universal phase shift. 

1.12 Nuclear spin state manipulation 

As mentioned in the previous sections, for a system of unpolarized nuclei 
and an effective hyperfine interaction strength A, the dephasing time in a weak 
magnetic field is T| ~ l/o" ~ y/N /A, where a is the width of the distribution 
of the nuclear field hz- This decay finds its origin in the ensemble average 
over the field distribution. In order to prolong the electron spin coherence, 
narrowing of the nuclear field distribution was proposed in [28] as an alterna- 
tive to the strategy of polarizing the nuclear spins [22], that would require a 
polarization close to 100% to be efficient, which is currently not available [28]. 
Few methods for nuclear spin state narrowing have been studied, in Ref. [31] 
the narrowing is due to gate-controlled Rabi oscillations in a double quantum 
dot in which the exchange interaction oscillates, in Ref. [37] a scheme based on 
quantum phase estimation is envisioned for a single undriven spin in a single 
quantum dot, and in Ref. [32] the narrowing is achieved by optical preparation. 

1.12.1 Nuclear state narrowing by qubit state measurement. Here we 
discuss a nuclear state narrowing technique that has been proposed in [31]. 
Consider for definiteness the ESR Hamiltonian Eq. (j56p . The effective Zeeman 
splitting is given by 6 + /i", where b = g* hbBz and /i" is an eigenvalue of hz- 
The idea behind state narrowing is that the ESR driving give rise to the 
resonance condition h + — uj = 0, such that the evolution of the electron 
spin depends on the nuclear spin state and thus a determination of the electron 
spin evolution results in a determination of the nuclear spin state. 

The eigenvalues of the nuclear field, as already mentioned in the previous 
sections, are Gaussian distributed in equilibrium. The diagonal elements of 
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the nuclear spin density matrix are pi{h^,t = 0) = = exp(— (ft," — 

(^2))^/2<7^)/\/27ra", with mean {hz) and variance a. Therefore, initiahzing the 
electron spin in the state | j) at time t = 0, the probability to find the electron 
spin in the state | i) is given by 

P^{t) = j dKPi{K,Q)Pl{t), (63) 

where -P"(i) is the probability to find the electron spin in the state | |), for a 
given an eigenvalue of the nuclear field hz , 



pnt) = \{K\^{l\u^^''it)\'i)^\K)[ 

1 6? 



-'ac 



26lc + 4^^ I 



1-cos Q^/6i;t4^) 



(64) 



If at time t = t^we perform a measurement of the electron spin and find | j) , 
the diagonal element of the nuclear spin density matrix will change according 
to 

pj{h:,0) ^ p^j''^\h:,tm) = Pi{K,0)^^. (65) 

-a (Cm j 

In the case where a measurement is performed with a low time resolution At, 
At ^ 1/6, such that it gives the time averaged value, the probability turns 
out to be = \imT^oc{l/T) Jq dtPl'it) = blJ2{bl^ +4(5^). Therefore, a 
measurement on the electron spin with outcome | |) results in a multiplication 
of the nuclear spin density matrix by a Lorentzian, with width 6ac) centered 
around the that satisfied the condition b + — u = 0. The nuclear spin 
distribution, thus, undergoes a narrowing, resulting in an enhancement of the 
electron spin coherence, if b^c < cr. In the case that the measurement outcome 
is I t) the diagonal element of the nuclear spin density matrix will change 
according to 

pi{K,0) - pf'^\K,tm) = Pi(K'^) :^_p^lt2y (66) 

resulting in a reduced probability for the nuclear field to have a value that 
matches the resonance condition b + — u> = 0. 

This procedure can be iterated many times before changes due to the slow 
internal dynamics start to affect the nuclear spin state. Many measurement 
of the electron spin are possible within this time, with re-initialization of the 
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electron spin state between the measurements. Assuming that M cycles can 
be performed with a static nuclear field, we have 

where is the number of measurement outcomes | j). If the outcome is 
I 1) the narrowing has been achieved, otherwise, it is necessary to wait for a 
re-equilibration of the nuclear system before the next measurement. 



1.12.2 Optical preparation of nuclear spins. Here we discuss the case of 
optical nuclear spin preparation that makes use of spin-flip two-photon Ra- 
man resonance in a driven three-level system (TLS) [32]. The lowest electronic 
states in GaAs quantum dots that are optically active under (T+ circularly po- 
larized excitation are the ground state of a single localized conduction-band 
{Ec) electron, in which a Zeeman field splits the up and down spin states, 
and the negatively charged exciton (trion) \X), given by two electrons with 
antiparallel spin plus one valence band heavy hole (hh) with angular momen- 
tum Jz' = +3/2, as schematically shown in Fig. [H The J = 3/2 subspace in 
the valence band split up into heavy and light holes (hh and Ih) along the di- 
rection z' of strong quantum dot confinement, that is in general different from 
the z-axis in the conduction band, chosen to be the direction of the magnetic 
field B. The two circularly polarized lasers stimulate the transition between 
I t) and \X) at frequency ojp = ujx — a;| — Ai and the transition between | [) 
and \X) at frequency = ~ ~ ^2j while the trion J^/ = —3/2 is not 
excited. 

The narrowing of the nuclear field distribution is based on light scattering 
in the TLS, where two long-lived spin states are resonantly coupled to the 
excited trion state \X) that decays spontaneously. For the two-photon reso- 
nance condition 5 = Ai — A2 = 0, where 5 is the detuning of the difference of 
the frequency of the two lasers ojc — ujp from the Zeeman splitting ojz of the 
two spin states, the system is in a superposition of the two spin states with a 
vanishing excited state component, and the system is driven to a dark state 
with no photon emission. In presence of a nuclear spin field, the resonance 
moves to 6 = 6hz, where 5hz is the deviation of the Overhauser field from 
its mean. Monitoring the photon emission constitutes a continuous weak mea- 
surement of the Overhauser field hz. The absence of photon emission in the 
limit t — > 00, corresponding to the strong measurement limit, would project 
the nuclear state onto \5hz = 0), with width a = 0, therefore letting the de- 
phasing time to diverge, T| ~ l/o" ^ 00. A continuous weak measurement of 
the Overhauser field, supported by an adaptive adjustment of the lasers fre- 
quencies every time a photon is detected, leads to a narrowing of the nuclear 
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Figure 4. Three-level system. The states 1(2) are spin-up (-down) conduction-band (Ec) electron, 
with splitting given by g/iB Btot + <5/iz, where Shz is the z component of the nuclear field 
fluctuations. State 3 is a trion with J^i =3/2 (With permission from the Authors). 

field distribution, and an enhancement of tfie pliase cofierence of the electron 
spin. 

The relevant effective Hamiltonian of the TLS in the rotating wave approx- 
imation is block diagonal, with blocks labeled by the eigenvalues of the 
field 5hz 

^ /dh'; + 6 np\ 
Wfc = -- Qc , (68) 

\ np -A J 

where A = Ai + A2. The combined system consisting of the TLS and the 
nuclear spins evolves in time according to a generalized master equation p = 
Cp, where C is the Liouvillian operator defined as 

p = Cp=l[H,p]+Wp, (69) 
in 

where H is the Hamiltonian of the system, and W = ^ ^Xai'^^^aXP'^Xa — 

(Txxp- P<yxx)/'i + Y.i3=i,x^pi'^^PPP'^Pp-'^PPP- P'^pp)/'^- Here the rate Txa 
describes the radiative decay of |X) into a = | j), | |), while 7/3 is the pure 
dephasing rate of state /? = | |), \X) with respect to | j). 

Taking as initial state a product of arbitrary density matrices xo for the 
TLS and = X^fcfc' ^kk'\Sh^){hz \ for the nuclear field, the stationary solution 
is an entangled state p. 

In order to describe the state of the system conditional on a measurement 
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record, a conditional density matrix is used. The a posteriori distribution v^^ 
is found to be concentrated around the two photon resonance. The stationary 
emission rate is 

Tern = TiSpit) = TY,T^kk{X\pkk\X), (70) 

k 

where T = Fxt + Txj, and S is the cohapse operator, describing spontaneous 
emission of the state \X) into | j) and | j) at rates Tx'\ and Txi- the update 
rule for u upon photon emission is 

/ _ ^kk{X\pkk\X) , . 

The population in the Overhauser field 5hz corresponding to the two-photon 
resonance 5hz = 5 is depleted by the photon emission. The electron spin co- 
herence is quantified by the time dependence of the transverse electron spin 
component, which in turn is given by the Fourier transform of the nuclear 
field distribution, (5+(t)) = (?i/2) z^/cfe exp(iM/i^). The repeated observa- 
tion of the quantum dot photon emission and consequent adaption of the laser 
frequencies after each photon detection leads to a narrowing in the nuclear 
distribution and consequent enhancement of the electron spin coherence time. 

1.12.3 Exponential decay in narrowed nuclear state. We have seen that 
the nuclear spin bath induces a non-Markovian dynamics of the electron spin, 
with super-exponential or power-law decay of the correlation functions. On the 
other hand it has been argued that a narrowing of the nuclear spin distribution 
is expected to prolong the electron spin coherence. In [38] it is shown that, 
in case of a large Zeeman splitting and for a particular narrowed nuclear 
spin state, a Markovian dynamics can arise from virtual flip-flops between the 
electron spin and the nuclear spin system, with simple exponential decay. 

The Hamiltonian describing the interaction of the electron spin with the 
nuclear system in a large magnetic field is given by Eq. ()53p . The energy non- 
conserving term V can be eliminated by means of a Schrieffer- Wolff transfor- 
mation, Ti = e^7ie~^ k, H = Tio + ^[S, V], where S = CqV, and Cq is the 
unperturbed Liouvillian, defined by CqO = [Ho, O]. The effective Hamiltonian 
H is given by 

H={lo + X)S' + D. (72) 

The operators lo, D and X are nuclear spin operators and the first two are 
diagonal in a product-state basis of /^-eigenstates, whereas X is purely off- 
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diagonal and produces correlations between nuclear spins. Corrections of the 
order of ~ /Nb in the diagonal terms of H are neglected, whereas the term 
of this size in X are retained. This assumption is valid as long as the bath 
correlation time Tc is much shorter than the time scale after which the diagonal 
corrections become relevant for 6 3> ^, where a Born-Markov approximation 
is valid: Tc ~ N/A < Nb/A^. As a result w = 6 + 

The electron and nuclear states are assumed to be initially unentangled 
and the nuclear system is prepared in a narrowed state, uj\n) = a;„|n). For 
this initial conditions, the dynamics of the transverse electron spin component 
is described by a GME, and can be written in a rotating frame defined by 
Xt = exp[— i(t(;„ + Aci;)t](S'+)t, where Aoj is a frequency shift sclf-consistently 
defined by Aw = -Re dtT.{t), with t{t) = exp[-z(a;„ + Aa;)i]S(t), through 
the memory kernel of the GME. The equation of motion for xt is given 
by 



Xt 



\ f drt{T)xt-r. (73) 



If decays to zero sufficiently fast on the time scale ^72, where T2 is in 
turn the decay time of xt, it is possible to approximate xt=T ~ xt and extend 
the upper limit of the integral to infinity, t — 00, obtaining a Markovian 
dynamics 

= exp(=t/T2)xo + e(t), — = -Im / dtT.{t), (74) 

-t2 Jo 

where e(t) gives a small non-Markovian correction that cab be bounded pre- 
cisely if S(t) is known. 

For a homo nuclear system, by expanding E(t) in the perturbation V = XS^ 
and retaining only leading orders in the Born approximation in the small 
parameter A/u^ the decoherence time T2 can be cast in the compact form 



T2 



coo 

Re 



poo 

/ dte-'^^\X{t)X), X{t)=e-''^*Xe''^\ (75) 
Jo 



where the average stands for an expectation value taken with respect to the 
initial nuclear state. Though the compact form resembles the standard result 
for pure dephasing valid in a weak coupling expansion, here there is no such 
weak coupling expansion. The decoherence rate I/T2 depends on the correlator 

C{t) = {X[t)X). For an isotropic electron wave function of the form '4){r) = 
'ip(0)e^^^^^"^''^'^ containing A'^ ^ 1 nuclei within a radius tq in d dimension, the 
asymptotic dependence of C{t) at long times is C{t) oc l/t^'^/'i, for t N/A 
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and d/q < 2. For d/q < 1/2, I/T2 given by Eq. (j75p diverges and no Markov 
approximation is valid within the Born approximation. On the other hand, for 
a 2D dot with a Gaussian electron wave function and for unpolarized nuclear 
system, Eq. (j75|) gives the simple result 



The condition for the validity of the the Markov approximation, T2 > Tc ^ 
N/A is satisfied whenever A/b < 1, which correspond to the range of validity of 
the Born approximation. Remarkably, from last equation it follows that I/T2 
strongly depends on the magnitude of the nuclear spin, I/T2 oc I^. Therefore, 
systems with large nuclear spin, such as In (/in = 9/2), will show faster decay. 

With these last results on exponential decay in a spin bath, we conclude 
the part on electron spin decoherence induced by the nuclear spin system and 
focus on phonon-mediated relaxation of the electron spin in quantum dots. 

1.13 Phonon-induced relaxation in quantum dots 

Electron spin relaxation in quantum dots takes place via transitions between 
spin states, with consequent energy dissipation in the environment. In a quan- 
tum dot the dissipative environment is represented mainly by the phonons 
in the surrounding crystal. Therefore, in order to fully understand relaxation 
and decoherence mechanisms that occur in quantum dots, it is important to 
understand the manner in which the electron spin interacts with phonons. 

Spin-orbit interaction creates an admixture of orbital and spin degrees of 
freedom of the electron, and represents an effective coupling mechanism that 
mediates the spin-phonon interaction, and that, ultimately, is responsible for 
relaxation of the electron spin. Phonons can produce electric field fluctuations 
that can lead to spin relaxation of eigenstates of the spin-orbit Hamiltonian. 
Two kinds of electron-phonon interactions are taken into account, that arise, 
respectively, from an inhomogeneous deformation of the crystal potential, re- 
sulting in an alteration of the band-gap, and a homogeneous strain due to 
piezo-electric effect, the former taking place in all semiconductors, the latter 
only in crystals without structure inversion symmetry such as GaAs. 

1.14 Introduction: Spin-orbit interaction 

An electron that moves in an electric field experiences an effective magnetic 
field in its rest frame which interacts with the spin of the electron. The internal 
magnetic field depends on the orbital the electron occupies and therefore spin 




(76) 
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and orbit are coupled. This well known effect comes directly from the relativis- 
tic Dirac theory of point particles and it goes under the name of spin-orbit 
(SO) interaction. The SO Hamiltonian has the general form [39] 

^so = — • {cT X VF), (77) 

where ttiq is the free electron mass, c is the speed of light, a = (cJx, fy, Cz) 
is the Pauli matrix vector, and V is the electric potential. In presence of an 
external magnetic field B = V x A, the canonical momentum p is replaced by 
the kinetic momentum P = p + eA, A being the vector potential. 

In semiconductors like Si or Ge the crystal lattice has spatial inversion sym- 
metry. For such materials, states of a given momentum k are 4- fold degenerate 
at i? = 0. In fact due to time reversal symmetry, ek,-f = e-k,j, holds, and from 
the inversion symmetry one has ek,o- = e-k.o-, such that ek,! = ek,j, = holds. 

The double degeneracy can be broken either via the application of an exter- 
nal magnetic field, which breaks the time reversal symmetry, or via the brake 
of spatial inversion symmetry. This is indeed what happens in crystals that 
exhibit bulk inversion asymmetry (BIA), such as the zincblende structure of 
GaAs. This effect is know as Dresselhaus spin-orbit interaction [40,41]. The 
Hamiltonian for 2D systems results from the 3D bulk Hamiltonian [42] after 
integration over the growth direction z along [001] 

Hd OC [-(J.j:p.j;{pl) + (JyPy{pl) + a^Pxpl - CTyPyPl] (78) 

where x and y point along the crystallographic directions [100] and [010]. Due 
to the strong confinement along z, the terms cubic in momentum components 
appearing in the Hamiltonian are usually much smaller than the linear ones, 
and they are usually neglected. Retaining only the linear term 

'Hd = (3{-(^xPx + CTyPy) (79) 

where (3 depends on material properties and on {pD- The spin dynamics re- 
sulting from the Dresselhaus Hamiltonian is well understood in the case of 
circular orbit, in which the spin rotates in the opposite direction with respect 
to the orbit, as shown in Fig. [3 

In heterostructures like GaAs/AlGaAs, an asymmetric confining potential 
additionally breaks the inversion symmetry, giving rise a further spin-orbit 
interaction due to structural inversion asymmetry (SIA), known as Bychkov- 
Rashba term [43,44]. For a confining electric field along the z direction, the 
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Figure 5. Schematic representation the apparent momentum dependent field f2(p) in the 
spin-orbit Hamiltonian, fJ(p) ■ <t, for the Dresselhaus and Rashba spin-orbit interactions. 



Bychkov-Rashba Hamiltonian oc (E x p) • <t is 



where a depends on the confining potential and on material properties. The 
spin dynamics resulting from the Bychkov-Rashba Hamiltonian can also be 
well understood in the case of circular orbit, in which the spin rotates along 
in the same direction as the orbit, being the spin always antiparallel to the 
direction of motion, as explained in Fig. [5j 



1.15 Electron spin relaxation and decoherence 

Due to the fact that the most promising semiconductor devices that make use 
of the electron spin as the quantum two level system (qubit) are realized on the 
basis of 2D electron gases in GaAs heterostructures, the following discussion 
concentrates on spin flip mechanisms that are relevant for GaAs. 

Spin relaxation of localized electrons in quantum dots shows remarkable dif- 
ferences from the case of delocalized bulk electrons. The most effective mech- 
anisms in bulk 2D are related to the broken inversion symmetry, either the 
BIA or SIA case, which give rise to a strong spin-orbit splitting in the elec- 
tron spectrum, ultimately responsible for spin flip. Besides, the piezoelectric 
interaction arising in non-inversion symmetric crystals provides a strong cou- 
pling of electrons to the bosonic bath of acoustic phonons. The interplay of 
these mechanisms results in an efficient spin-lattice relaxation for bulk carriers 
in III-V type semiconductors and heterostructures. The strong localization of 
electrons in quantum dots leads to suppression of spin-flip rate. The phonon- 
assisted spin-flip mechanisms in semiconductor quantum dots have been stud- 
ied in [18,19]. 




(80) 
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1.15.1 Electron spin relaxation in quantum dots. In the case of strong 
confinement in the z direction, corresponding to the [100] crystallographic 
axis, for a lateral dot size much larger than the degree of vertical confinement, 
the Hamiltonian derived from the Kane model [45] for 2D electrons in the 
conduction band in the presence of an external magnetic field B is [18, 19] 

W = 1^ + U{r) + t/ph(r, t) + ^gfiBCT • B + + ^SO' (81) 
'^SO = Pi-CTxPx + CTyPy), T-^so = aicr^Py - CFyPx). (82) 

Here p = —ikV + (e/c)A is the kinetic momentum, m the effective mass, g 
the effective electron g(-factor (in GaAs g = —0.44), and a the Pauli matrix 
vector. The axes x, y, and z coincide with the main crystallographic ones, 
with z perpendicular to the 2D plane. The first two terms of the Hamiltonian 
describe the quantum dot with confining potential U{r), that is typically cho- 
sen parabolic. The third term describes the spin-independent interaction with 
acoustic phonon. The fourth term is the Zeeman Hamiltonian. TCso describes 
the spin-orbit effects. Ti-gQ is the Dresselhaus term, due to BIA, and Ti-gQ is 
the Rashba term, due to SIA. For GaAs heterostructures /3 ~ 10^ cm/s. 

The Hamiltonian of Eq. (jSip should contain also term describing "direct" 
interaction between spin and phonons, such as that due to an inhomogeneous 
deformation of the lattice, and a term describing the spin-phonon coupling in 
presence of a magnetic field due to a lattice-deformation-dependent admixture 
of valence-band and conduction-band states. Their contribution on spin relax- 
ation rates turns out to be negligible with respect to the dominant admixture 
mechanism contribution ascribable to the Dresselhaus and Rashba spin-orbit 
interaction. See Ref. [18,19,46] for a discussion of the direct spin-phonon cou- 
pling contribution. 

The phonon-induced rate for the transition between |^'n) and is given 
by Fermi's golden rule 

r = Yl2\(^n\n,H\^i)\'D{ez). (83) 

n 

Here D{E) is the phonon density of states at the Zeeman energy splitting e^. 
From experimental results, the relevant acoustic phonons can be treated as 
bulk-like phonons, showing a linear dispersion relation in the relevant energy 
range, for which the density of states increases quadratically with energy [20]. 
The states and \^n) are the effective spin states, containing more than 
one orbital and both the spin up and down states. This admixture of spin and 
orbit comes out in taking into account the spin-orbit interaction due to BIA 
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Figure 6. Schematic representation of the spin-flip process associated with the emission of a 
phonon of energy tz and momentum qj, where j is the branch index, at rate Fi. 



and SIA as a perturbation. Due to the localization of stationary states in a 
quantum dot, it follows that the spin-orbit interaction does not directly couple 
Zeeman-split sublevels in the same quantum dot orbital. It follows that within 
first order perturbation theory in the spin-orbit Hamiltonian, the effective 
single electron quantum dot states are 

i<>-i"T)+E ^„il:°5L^ Ki>, (84) 



where {Tisohn'ri ~ ('^' i \'^so\f^ T) and{|n)} are the unperturbed quantum dot 
orbital states. Due to the anisotropy of BIA and SIA spin-orbit interaction, the 
admixture of spin and orbit degrees of freedom turns out to be anisotropic [47] . 

For spin-flip transitions involving a small energy transfer, the dominant con- 
tribution comes from piezo-electric phonons. The electric field associated with 
a single phonon scales like 1 / -^g for piezo-phonons and like for deforma- 
tion potential phonons, q being the phonon wave number. This is due to the 
fact that piezo-phonons come from a homogeneous lattice strain, in which long 
wavelengths play a major role. Vice versa, a local deformation would involve 
short wavelengths, and so higher energies. On the other hand, wavelengths 
much longer than the dot size give rise to a global shift of the entire dot po- 
tential, therefore the effective phonon wavelengths are those comparable with 
the dot size, as seen in [48]. 

The electron-phonon coupling for piezo-electric phonons has the form 



U^l^{r, t) oc exp(fq • r - iuJqat)Aciab\a + ^-C-, 



(86) 
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where qa are the wavelength and the branch of the phononic modes, A^^ is 
effective anisotropic piezo-electric modulus of wave qa. The matrix elements 
of the phonon Hamiltonian between the Zeeman split sublevels of orbital level 
n, that describe the spin-flip process with emission of phonon qa, are given 
at first order in the spin-orbit interaction by 



En — Ek — giisB En — Ek+ gfisB 



(87) 



As a consequence of Kramer's theorem, in case of no external magnetic field, 
Eq. (|87|) is zero. Considering only Dresselhaus spin-orbit interaction, for small 
Zeeman splitting, g^sB <^ \/ms'^huiQ, it is possible to obtain an effective spin- 
flip Hamiltonian which acts on the subspace of Zeeman sublevels of orbital level 
n, [19], where a phonon induced electric held arise as a gradient of Uph. For 
a parabolic dot confinement potential, and for the particular case of circular 
dot with level spacing ujq, the spin- flip rate for the transition between the 
Zeeman sublevels of the dot ground state, associated with the emission of a 
piezo-phonon as depicted in Fig. [6l is [19] 

^ _ {gfiBB)\ , o,^,,, , _ 2 (e/ii4)2/32 / 1 



A,(l + cos^(^)), A, = — + 



where (3 is the strength of the Dresselhaus spin-orbit interaction, i? is the angle 
between the direction of conflnement in the quantum dot z, and the direction of 
the applied magnetic held z', and Ap is the strength of the effective spin-piezo- 
phonon coupling. For given longitudinal and transverse sound speed si and st, 
crystal mass density /j, and modulus of the piezo-tensor e/ii4 (e/114 = 1.3 x 10^ 
eV/cm for GaAs), it ranges from 7 x 10~^ to ~ 6 x 10~^, depending on (3. 
Eq. ([88|) shows a strong dependence on the magnetic held and the lateral dot 
conflnement energy (Jq. For Itloq = 10 K and magnetic field i? = 1 T, Fi ~ 
1.5 X 10^s~^. These theoretical expectations, in particular the i?-dependence, 
have been confirmed in experiments [49-51], and long spin relaxation time, up 
to Is, have been measured [21]. 

The effect of the Rashba spin-orbit interaction has so far not been taken 
into account. As the Dresselhaus term, it contributes to the admixing of spin 
and orbital states, and therefore to relaxation due to phonon scattering. The 
effect of the interplay of these two terms can show up in a strong difference in 
the their associated relaxation rates [21,52-54]. For a quantum dot in external 
magnetic field, the first and second lowest levels show a crossing behavior as 
a function of the applied magnetic field, the ground state being not affected. 
In a perturbative treatment of the spin-orbit interaction [52], Dresselhaus and 
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Figure 7. Schematic representation of the relaxation process (left side), due to magnetic field 
fluctuations 5B^{t) orthogonal to the applied B field, and of the dephasing process (right side), 
due to magnetic field fiuctuations (i) parallel to the applied B field. 



Rashba terms show a qualitative difference, in which the latter couples the 
crossing levels, giving rise to an anticrossing of the levels at the point of ac- 
cidental degeneracy. For magnetic fields much smaller than the crossing level 
value, the these two levels have a well defined spin orientation, i.e. low degree 
of admixture. In the region of anticrossing the admixture leads to complete 
superposition of the up and down states, and eventually to a reversed situation 
in the limit of magnetic field much larger than the crossing value, in which 
the two levels have again well defined spin, but reversed. Therefore, sweeping 
the magnetic field over the crossing region leads to spin-flip. In particular at 
the avoided crossing point, the strong admixture between spin states lead to 
a cusp-like behavior of the relaxation rate as a function of the magnetic field, 
at the anticrossing point. 



1.15.2 Phonon-induced electron spin decoherence. In a Markovian dy- 
namics the decoherence time T2 is limited by both spin-flip and dephasing pro- 
cesses, though its upper bound is T2 < 2Ti. A systematic analysis of phonon- 
induced spin decay is carried out in [55]. Also there both the Rashba and Dres- 
selhaus spin-orbit interactions are treated perturbatively. The deformation 
potential phonons are also considered. For a spin-orbit length A50 = fi/fn* (5 
much larger than the electron orbit size A, the contribution to the spin-phonon 
coupling in the Hamiltonian linear in X/\so oc a, /3 is due only to a flnite Zee- 
man splitting. For B field in the range m*P'^ ^ gfJ-BB <^ ?k<Jo, the effective 
Hamiltonian is 

nes = ^gt^B[B + 5Bit)]-(T, 6B{t) = 2Bxn{t), (89) 

where n{t) = (-01 [(L^^^), [/ph(t)]|^/^), lip) is the electron orbital wave function, 
Lrf is the dot Liouvillian, L^A = [fid, A]. The vector ^ lies in the 2D dot plane 
and depends on a, /?, and m* . The most important consequence of Eq. (j89p is 
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that within first order in the spin-orbit couphng parameter there can only be 
transverse fluctuations of the effective field. 

In case of many uncorrelated scattering events, the expectation values (S) 
of the spin obeys the Bloch equation [6] 

(S) =a> X (S) -r(S) + T, (90) 

where u = u'B/B, uo = gfisB/h. In the Born-Markov approximation, for a 
generic 5B such that ((5B) = 0, the tensor T can be expressed in terms of the 
spectral function 

Mw) = ^ dt{6Bm^Bj{t))e-'^\ (91) 

and it is diagonal in a frame {X, Y, Z), with Z oriented along B. The symmetric 
part of r is responsible for the decay of the spin components, and it can be 
expressed just as function of jfj{w) = Re[Jij(t(;) ± Jij{—w)\. T can be split 

in two contributions, T = + T'^, where contains the spectral function 
Jij(w) at the Zeeman frequency oj, and describes spin decay due to emission 
or absorption of a phonon, whereas is due to elastic scattering of spin. The 
relaxation time Ti is completely determined by r**, while the decoherence time 
T2 is affected by both T^ and F*^, the latter describing pure dephasing. In the 
general case, solution of Eq. (190]) yields 

^^■.= Tzz = T^zz. y^--=\{Txx+Tyy). (92) 

In many cases the contribution of V' to spin decoherence is negligible, the 
decoherence rate being determined entirely by V^. However, it turns out that, 
at first order perturbation theory in the spin-orbit interaction, no dephasing 
takes place [55]. Due to the transverse nature of the fluctuations in the mag- 
netic field in the effective Hamiltonian, the tensor F*^ is identically zero, F'^ = 
(Fig. [7] illustrates an intuitive picture of the effect of the longitudinal and the 
transverse fluctuations). As a result 

^ = ^ = 4xM + '^^yM- (93) 

Contributions to the decoherence time T2 due to pure dephasing arise when 
two-phonon processes are taken into account in the next order in the electron- 
phonon interaction [56]. Therefore, if only spin-orbit decay mechanisms are 
taken into account, the decoherence time T2 for the decay of the transverse 
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component of an electron spin in GaAs quantum dots is T2 = 2Ti . 
1.16 Spin-orbit interaction for heavy holes 

The electron spin in GaAs quantum dots has shown to have a long relaxation 
time, due to inefficient phonon-induced relaxation mechanisms. On the other 
hand, the decoherence time is mainly dominated by hyper fine induced decay, 
due to the fact that the decay of the longitudinal electron spin component 
can be strongly suppressed by the application of an external magnetic field. 
In order to circumvent this problem, the use of hole spins as carriers has 
been recently proposed. The valence band in III-V semiconductors has a p 
symmetry, for which the electron has zero probability to be found on the 
position of the nucleus. According to Eq. ()5ip . it follows that the hyperfine 
interaction between holes and nuclei is strongly suppressed with respect to 
that of nuclei and conduction band electrons. However, the hole spin relaxation 
time turns out to be much smaller than that of the electrons by several order 
of magnitude. The reason for this is due to the fact that, beside the spin- 
orbit coupling due to bulk inversion asymmetry and the structural inversion 
asymmetry, there is a strong spin-orbit coupling between the heavy-holes (HH) 
and the light-holes (LH) sub-bands [57]. Investigations of hole spin relaxation 
in quantum dots, exclusively due to spin-orbit coupling of LH and HH sub- 
bands, give estimates for the relaxation time much shorter than the case for 
electron spin [58,59]. 

In Ref . [64] , HH spin relaxation is analyzed in presence of Rashba and Dres- 
selhaus spin-orbit coupling, as well as spin-orbit between HH and LH. From 
the two-band Kane model, the Hamiltonian for the valence band of III-V semi- 
conductors is given by 

Wbuik = WLK + r/J-n + Wz, (94) 

where Wlk is the Luttinger-Kohn Hamiltonian [60], rj oc {Eg + Aso)/Aso, Ago 
is the split-off gap energy, and Eg is the band gap energy. J = {Jx,Jy,Jz) 
are the 4x4 matrices corresponding to spin 3/2, Jl^ = Pz{Px ~ Py)i and 
Qx, are obtained by cyclic permutations. The last term in Eq. (j94p Tiz = 
—2k,ijlb^ ■ J — 2qfiB^ ■ is the Zeeman term for the valence band [61], with 
K and q Luttinger parameters [61], and J' = {J^, Jy , Jz)- 

In case of structure inversion asymmetry along the growth direction, due to 
an asymmetric confinement, there is an additional contribution to the spin- 
orbit interaction, the Bychkov-Rashba term. For the two-band Kane model it 
is given by [62,63] aji(P x E) • J, where E is the effective electric field along the 
growth direction, and a/j is the Bychkov-Rashba spin-orbit coupling constant. 
We consider a two-dimensional system grown along the [001]-direction. Be- 
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cause of confinement, the valence band splits into a heavy- hole subband, with 
Jz = ±3/2, and a light-hole subband, with = ±1/2 [57,64], where z is the 
growth direction. In case of large splitting A between HH and LH, the proper- 
ties of the two subband s can be described separately, the Jz = ±3/2 subspace 
for HHs and the Jz = ±1/2 subspace for LHs, using only the 2x2 submatrices. 
The HHs submatrices have the properties Jx = Jy = 0, and Jz — ^ctz [65]. For 
low temperatures only the HH subband is significantly occupied. Considering 
only HHs, starting from the bulk Hamiltonian Eq. (j94p with the addition of 
the Bychkov-Rashba term, at the lowest order in perturbation theory [66], it 
is possible to derive an effective Hamiltonian for a quantum dot with lateral 
confinement potential U{x,y) 

n = \{Pl + P^) + U{x, y) ± - ^gzzfisBzaz, (95) 

where m is the effective HH mass, gzz is the component of the g factor tensor 
along the growth direction, and the effect of an in-plane component of the 
magnetic field can be neglected due to strong anisotropy in the HH g factor, 
g\\ ^ gzz [65], as well as the orbital effect of the in-plane magnetic field, as long 
as S|| <C ch/eh'^, h being the height of the quantum dot. P = p± (|e|/c)A(r), 
with A(r) = {-yBz/2,xBz/2,yBx - xBy), and 

nf§ = iaa+P^ + (3P-P+P-(T+ + 7S_PV^ _^ j^.c. (96) 

The first two terms in the spin-orbit interaction for heavy holes consist in the 
Rashba and Dresselhaus contribution, respectively, while the last term (7) de- 
scribes the combination of two effects: the orbital coupling via non-diagonal 
elements in the Luttinger-Kohn Hamiltonian (oc P±), taken into account per- 
turbatively, and magnetic coupling via non-diagonal elements in the Zeeman 
term, (oc B±) [61]. This new spin-orbit term is unique for heavy holes [67]. In 
Eq. ill) Q - 37oaK(^^)/2moA, (5 = -37or/(p2)/2moA, 7 = 37oK/XB/moA, 
a± = {cjx ± Cj;)/2, P± = Px^ Py, and B± = B^ i: By, mo is the free electron 
mass, 7o is the Luttinger parameter [61], {Ez) is the average electric field, and 
A is the splitting between heavy and light hole subbands, A tx where h 
is the quantum-dot height. For a quantum dot with characteristic lateral size 
I, the ratio (Wso)/(^so) (V^)^- Therefore for flat quantum dots, l/h » 1, 
the spin-orbit coupling for heavy holes can be weaker than that for conduction 
electrons [52,55]. This observation has also been confirmed experimentally [68], 
where the spin relaxation rate for heavy holes has shown to be comparable to 
that of electrons. 

For vanishing spin-orbit interaction, the spectrum of the Hamiltonian Eq. 
(195p for a parabolic lateral confinement can be found through a canonical 
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transformation [69], and it is the Fock-Darwin spectrum spht by the Zeeman 
term [70,71]. In the framework of perturbation theory [52], it can be seen that 
the corrections to the spectrum due to arise only at second order, and 
the spin-orbit interaction influences the wave functions more strongly than the 
energy levels. couples the two lowest states |0, ±3/2) to the states with 

opposite spin orientations and different orbital momenta |/,=f3/2). The dif- 
ferent spin-orbit interactions appearing in Eq. (j96|) differ by symmetry in the 
momentum space [52,72], and thus produce a mixing of spin-up and spin-down 
states, with resulting avoided crossings between energy levels. We mention here 
that the levels cross only if g^z > 0, therefore in case of GaAs (qzz > 0) quan- 
tum dots an anticrossing appears, with consequent peak of the relaxation rate 
as a function of the magnetic field, at the point where the crossing takes place, 
while for InAs {gzz < 0) quantum dot no crossing and no cuspic-like behavior 
of the relaxation rate appear. The spin-orbit mixing of the heavy-hole states 
provides a mechanism of transitions between the states |0, ±3/2) through emis- 
sion or absorption of an acoustic phonon, that ultimately represents the main 
source of relaxation and decoherence for heavy- holes [64]. 

Taking into account piezoelectric and deformation potential phonons, the 
potential of a phonon with mode qa is given by [19, 66] 



h 



2pSaqV 



FiQz 



a + 



-MzdT 



(97) 

where qy = {qx,Qy), o, and b are constants of the deformation potential, V 
the quantum dot volume, Sa the sound velocity, p the crystal 
effective piezoelectric modulus, dqc, the phonon polarization vector, F{qz) the 
form factor, which is determined by the spread of the electron wave function 
in the z direction. 



1.16.1 Spin decoherence and relaxation for heavy holes. For a single- 
particle quantum dot, in which an heavy hole can occupy one of the low- lying 
levels, some energy levels with same spin orientation cross, with increasing B, 
the upper Zeeman-split ground state level. Therefore we consider an n-level 
system, in which the first n — 1 levels have same spin orientation, while the 
n-level has opposite spin. In the context of Bloch-Redfield theory, the Bloch 
equations for the spin motion of a heavy hole in such a system are given, in 
the interaction picture, by 



{Sz) = {St - {Sz))/Ti - R{t), (98) 
{S.) = -{S.)/T2, {Sy) = -{Sy)/T2, (99) 
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where R{t) = WniPnn{t) + Y17=i ^niPuit), p{t) is the density matrix, Wij is 
the transition rate from state j to state i, is a constant that takes the value 
{Sz) in the thermodynamic equiUbrium R{t) = 0, and 

1 111 

i=l i=l 

where the pure dephasing term, which is due to fluctuations along the longi- 
tudinal z direction, is absent in the decoherence time because the spectral 
function is superohmic. The spin motion involves n — 1 states and therefore 
there are n — 1 transition rates. It can be shown, by solving the master equa- 
tion, that for low temperature hqa ^ ksT, R(t) ^ and phonon absorption is 
strongly suppressed. In this case only one relaxation rate contributes to relax- 
ation time Ti. In this limit the last sum in the expression for the decoherence 
rate can be neglected, and the decoherence time T2 saturates, T2 = 2Ti. The 
relaxation rates for the different spin-orbit interactions are [64] 



A^o^aWj-^ , ) , (101) 

±^oc(3'u;l(^ (102) 



jmI II \2oj++ujz 2u)- — uz 



where uj± = $7 it uJc/2, ^ = \J + uj1/A uj^ = \e\B/mc is the cyclotron 



frequency, ujz — SzzPbBz, and i^y — y + By. In contrast to the case of 

conduction electrons [55], no interference takes place for heavy holes, and the 
rates originating from different spin-orbit terms sum up, giving the total spin 
relaxation rate 1/Ti = 1/Tf ^ + 1/Tf + I/t]'. For the case of GaAs quantum 
dots the crossing between levels takes place at uJz = ^-i 2i^_, and 3a;_, and 
the strong spin mixing arising causes cuspic-like peaks in the relaxation rate 
as a function of the external field B. 



1.16.2 Electric dipole spin resonance for heavy holes. The possibility of 

coherent manipulating the spins is of great importance for spintronics and 
quantum computation. In case of conduction electron spin-based electronics, 
such control is obtained by the electron spin resonance (ESR). Through the 
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application of short resonant microwave pulses, arbitrary superpositions of 
spin-up and spin-down states can be created. Therefore ESR provides a nec- 
essary tool for single-qubit operations, an essential requirement for quantum 
computation. In Rabi oscillations and spin echo experiments [73], that are 
based on this technique, the ESR signal can be detected by measuring the 
absorption of radio- frequency (rf) power [74]. ESR methods involve magnetic- 
dipole transitions induced by oscillating magnetic fields. Besides, an alternative 
is provided by alternating electric fields, that give rise to electric-dipole spin 
resonance (EDSR). 

Considering the spin-orbit coupling as a perturbation, at first order the 
two states corresponding to the Zeeman-split ground state |ib) are given as a 
superposition of few unperturbed Fock-Darwin states and spin states, |n, i)\s), 
with n G N the principal quantum number, \£\ < n the azimuthal quantum 
number, and s = ±3/2, for a detailed see Ref. [67]. In the case of heavy holes 
it can be shown [67] that magnctic-dipolc transition (An = 0, Ai = 0, and 
As = ±1) are forbidden, while, because of spin-orbit coupling between states 
with different orbital momenta and opposite spin orientations, |0, 0, ±3/2) and 
|1, ±1, ^3/2), electric-dipole transitions (An = ±1, Ai = ±1, and As = 0) are 
most likely to occur. Heavy holes are thus affected by the oscillating electric 
field component, but not by the magnetic one. EDSR for heavy hole appears to 
be an essential tool for the control of spin dynamics and for the determination 
of important parameters, as the effective g factor, effective mass m, spin-orbit 
coupling constants, and spin relaxation and decoherence time. 

The Hamiltonian for the interaction of HHs with a circularly polarized 
electric field, that rotates with frequency w in the Xy-plane, E(t) = 
E {sin cot, — cos ojt,0), is given hy = {\e\E /mu!){Px cos ujt + Py sin ujt). The 
coupling between the states |±) is given by = dsoEe~^^\ where 

dso = (|e|//2c^)(/3+cu+ +/3ru;_), (104) 

is an effective dipole moment of a heavy hole and it depends on Dresselhaus 
spin-orbit coupling constant, perpendicular magnetic field B_\_, lateral quan- 
tum dot size, and frequency u of the rf electric field. For details on (if and I 
see Ref. [67]. 

The effective master equation for the density matrix yO„^, in the contest 
of Bloch-Redfield theory, takes the form of Bloch equations, with a rf field 
detuned from uz, (^rf = — Larmor frequency 2dsoE/h, spin relaxation 
time Ti = l/(IV+_ -|- Wnm being the transition rate from state m to 

state n, decoherence time T2 = 2Ti, and equilibrium value of Pz without rf 
field given by = (VF+_ - W^+Ti). 

The coupling energy between a heavy hole and an oscillating field is given 
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by 

{H^(t)) = -dso • E(t), (105) 

where dso = dso{ip h — iPA + P h; 0) is the dipole moment of a heavy hole. 

The rf power P = —d{7i^{t))/dt = —udsoEp- absorbed by a heavy hole spin 
system in a stationary state is given by [75] 

_ 2u;{dsoEfT2p^/h 

l + 6lTi + {2dsoE/hYT,T2- ^ ^ 

The dependence of P on perpendicular magnetic field B± and frequency of 
the oscillating field ui shows three resonances and one resonant dip. The first 
resonance corresponds to the Zeeman energy of the heavy hole = tujj/g±^^ , 
the second to the first anticrossing between the unperturbed £"0,-3/2 £1,3/2 
energy levels, loz — ^jJ-, the third resonance reflects the peak in the decoherence 
rate due to an applied in-plane magnetic field at the second anticrossing 
LOz = 2a;_. The resonant dip takes place at zero dipole moment. 

The study of the position of these resonance allows to determine g±, m, 
and LdQ, while the shape and height provide information about the spin-orbit 
interaction constants a, (3, and spin-orbit interaction strength due to in-plane 
magnetic field. Besides, it is possible to extract informations about the depen- 
dence of spin relaxation and decoherence times on B±. 



1.17 Superconducting qubits 

Besides spin qubits in semiconductor quantum dots, superconducting qubits 
represent a category of promising candidates for the implementation of ar- 
tificial two-level systems as qubits. The key ingredient in building supercon- 
ducting qubits is the strong nonlinearity of the current-voltage relation of a 
Josephson junction. The ability to isolate few charge states on a supercon- 
ducting island, together with the possibility to let them interact through the 
coherent tunneling of Cooper pairs through the junction, represent a promis- 
ing way to control a operate a purely quantum system (charge qubits). The 
fiux quantization together with the strong nonlinear potential, arising from 
the current-voltage relation, provide a way to isolate few current states and 
coherently superimpose them (flux qubit). 

Superconducting qubits can be included in a more general framework of 
quantum circuits, that are electrical circuits showing, in the low temperature 
regime, quantum behavior, including quantum fluctuations [76]. In this con- 
text, as LC-circuits provide electrical realizations of quantum harmonic oscil- 
lators, showing a linear current-voltage relation, Josephson junctions provide 
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a full anharmonic counterpart, showing a more rich spectrum, with groups of 
few energy levels well separated from higher bands of the spectrum. 

Several types of superconducting qubits based on Josephson junctions have 
been so far theoretically proposed and experimentally realized (for compre- 
hensive reviews see [77,78]). Apart from the particular design of each device, 
superconducting qubits can be classified according of the working regime of 
the Josephson elements that constitutes the circuit. Every Josephson junction 
is characterized by two features: i) a critical current Ic, that is the maximal su- 
percurrent that can flow through the junction; and ii) an effective capacitance 
that the two superconducting faces have to accumulate charge. Together the 
energy associated with the critical current Ej = Ic$o/2vr and the charging 
energy of the associated capacitance Eq = e^/2C are the two most impor- 
tant parameters that determine the qubit working regime. For Ec 3> Ej the 
charge degrees of freedom are well defined and the number of Cooper pairs in 
a superconducting island is a well defined quantum number. Qubits that work 
in this regime are called charge qubits [79-88]. To the contrary, for Ec <C Ej 
flux degrees of freedom have well defined values, and current states are well 
defined. Qubits that operate in this regime are called flux qubits [89-92]. Other 
realizations of superconducting qubits, for different values of the ratio Ej / Ec, 
and many kind of possible accessible parameter regimes have been explored. 
The so called phase qubit [93, 94] operates in the flux regime, but is com- 
pletely represented by the superconducting phase, and it has no magnetic flux 
or circulating current associated. The quantronium [95], consisting of a split 
Cooper pair box arranged in a loop containing an extra large junction for the 
read-out. 

Experimental observation of Rabi oscillations in driven quantum circuits 
have shown several periods of coherent oscillations, confirming, to some ex- 
tent, the validity of the two-level approximation and possibility of coherently 
superimpose the computational two states of the system. Nevertheless, the 
unavoidable coupling to a dissipative environment surrounding the circuit rep- 
resents a source of relaxation and decoherence that limit the performances of 
the qubit for quantum computation tasks. Therefore, for the implementation 
of superconducting circuits as quantum bits, it is necessary to understand the 
way the system interacts with the environmental degrees of freedom, and to 
reduce their effect, if possible. 

1.18 Circuit theory and system Hamiltonian 

A systematic approach to obtain the Lagrangian of a generic circuit containing 
many different lumped elements, as well as Josephson junctions, has been 
proposed in [96,97]. In this way, it is possible to construct the full classical 
Hamiltonian of the system, quantize it and study its quantized spectrum, in 
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the two-level approximation. 

1.18.1 Network graph theory and the equations of motion. By means 
of classical network theory, an electric circuit is represented by an oriented 
graph, consisting of nodes and branches. Each branch correspond to a sin- 
gle two-terminal element, such as resistor, capacitor, inductor, current source, 
voltage source, etc.. The branches are then divided into two groups, the tree, 
representing a set of branches of the graph connecting all nodes without con- 
taining any loop, and the chords, represented by all the rest of branches. This 
way every time a chord branch is added to the tree a loop is obtained. The 
grouping in chords and tree depends on the formalism adopted, that in turn 
is functional to the kind of circuit described, either a flux qubit or a charge 
qubit. All the the topological information of the circuit is contained in the fun- 
damental loop matrix F, that connects tree branches and loops (i.e. chords), 
such that the matrix elements Fxy can be 1, -1, 0, depending whether the 
tree branch X and the chord branch Y have the same orientation, different 
orientation in the loop, or do not belong to the same loop. 

The equations of motion are represented by Kirchhoff 's laws, and can be at 
once written as 

Fleh = -Itr, F^Vtr = Veh - ^ext- (107) 

Here ^ext takes into account the possibility of having time dependent applied 
external fluxes. The fluxes and charges of the circuit represent the canonical 
variables of system, and they can be formally defined for the generic element 
X as 

Ix {t) = Qx (t) , Vx (t) = *x (t) . (108) 

Prom the last equation and from the second Josephson relation, it is possi- 
ble to identify the formal flux associated to the Josephson junction as the 
superconducting phase difference (f across the junction, 

where the second formula represents the first Josephson relation. With current- 
voltage relations for the various types of other branches, it is possible to obtain 
the classical equations of motion for the superconducting phases 



C(p = -LjhiiKf - Mocp -Md*<fi- — (N*ext + Sis), (110) 
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where Lj^ = 2it1c/^o is a diagonal matrix for the Josephson inductances of 
the junctions, Mq is the matrix of hnear inductance, describing their energy 
and mutual interaction, N and S describing the inductive coupling of the 
phases (p with external fluxes and currents, respectively. Md(t) is a symmetric 
matrix containing all the dissipative dynamics of (p, [96]. 

1.18.2 Two-level approximation. Dissipative elements present in the cir- 
cuit are incompatible with a Hamiltonian description of the system, therefore 
for the moment we omit them. In order to derive the Lagrangian for the elec- 
tric circuit, a complete set of unconstrained flux and charge degrees of freedom 
has to be isolated, such that every assignment of values to those charges and 
fluxes represents a possible dynamical state of the system. The Hamiltonian 
of the circuit follows straightforwardly from the Lagrangian by means of a 
Legendre transformation, and can be formally written as 

W = ^(Q - CvVf C-\Q - CvV) + (^)' (111) 

OyrJ . 1 27r 

Ui'P) = -y] cos + -<p^Mo<p + — <p^(N*ext + sis), (112) 

where C is the capacitance matrix, collecting all the capacitivc elements of the 
circuit, and describing the effective charge energy of the system, Cy describes 
the coupling of the charges Q to externally applied voltages V, The number 
of Cooper pair, that accumulates on a junction capacitance, and the phase of 
the superconducting order parameter through the junction, for sufficiently low 
temperatures, become quantized and satisfy canonical commutation rules, 

= ih5ij, (113) 

where 2e is the charge of a Cooper pair, and $o/27r is flux quantum. There- 
fore, once a Hamiltonian is obtained from circuit theory, its quantization fol- 
lows straightforwardly. The energy of an isolated system is a conserved quan- 
tity, therefore strictly speaking the Hamiltonian should be time independent. 
However time-dependent circuit elements, such as alternating currents and 
voltages, can be included in the Hamiltonian description as time-dependent 
parameters. 

Care should be taken when dissipative elements such as resistors are present 
in the circuit. In this case a more general approach must be adopted, in which 
the system considered is coupled to a environmental bath, and the dynamics of 



27r 



ipi,2eNj 
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Figure 8. Scanning electron micrograph of the Delft flux qubit (small loop with three Josephson 
junctions) and attached SQUID (large loop) [92] (With permission from the Authors). 

the circuit under analysis is obtain as the dynamics of part of a larger isolated 
system, as discussed in section 11.51 

Once the Hamiltonian has been obtained and quantized it is possible to 
study the temperature regime, in which few low energy states are taken into 
consideration. A two-level approximation can be carried out by considering 
only the ground state and first excited state, and neglecting higher levels of 
the spectrum. The goodness of the two-level approximation is controlled by 
the ratio of the temperature and the energy gap between the first and second 
excited state, k^T j IS.\2 ^ 1. The Hamiltonian of the two-level system can be 
therefore expressed in the form of a pseudo-spin 1/2 

n = ^(7, + (114) 

where A denotes the tunnel coupling between the two qubit states |0) and |1), 
eigenstates of cj^, and e the bias, due to asymmetry. 

1.19 Decoherence in superconducting qubits 

In this section we choose to describe decoherence effects in only two realizations 
of superconducting qubits, namely the phase qubit [93,94] and the flux qubit 
[89,96,98]. An extensive treatment of decoherence in the quantronium circuit 
is carried out in [99]. 

1.20 The superconducting persistent current qubit: Delft qubit 

In the working regime Ej ^ Eq, three types of circuit designs have been pro- 
posed, the Delft flux qubit [89,92], the IBM flux qubit [96], and its gradiometer 
variety [98]. The phase qubit operates in the same regime, see Sec. 11.21] 
The flux qubit realized at Delft [92] consists of a superconducting loop in- 
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terrupted by three Josephson junctions (see Fig. [8]). The strong flux regime 
Ej ^ Ec allows flux quantization of the flux through the loop, v?i + ¥?2 + V^s = 
27rn. Therefore, only two of the three phases of the junctions play the role of 
dynamical variables. For sufficiently low temperatures, in the small loop in- 
ductance limit, the inductive degree of freedom associated with the loop is 
frozen, and the effective potential U{(p) is periodic and shows a double well 
shape. The charging energy of the system here plays the role of the kinetic 
energy, and the Hamiltonian is written as 



The lowest energy states are two flux states localized in the two well minima 
(fii^ and (fji, and they correspond to clockwise and counter clockwise circulating 
currents in the loop, \L) and \R), encoding the logical |0) and |1) states of 
the qubit. Tunneling through the potential barrier between the wells lifts the 
degeneracy between the two current states, giving rise to a splitting A = 
{L\7i\R) between the lowest states of the system, that become the symmetric 
and symmetric superpositions of the current states. An external bias flux can 
create asymmetry in the double well, e = {L\7i\L) — {R\7i\R). The qubit 
Hamiltonian written in the {\L), \ R)} basis takes the form 



1.20.1 Markovian dynamics due to dissipative circuitry. The regime of 
working of flux qubits, in which the charging energy is much smaller than the 
Josephson energy, Ec <^ Ej, makes the flux qubits substantially insensitive 
to background charge fluctuations. Still, however, other mechanisms can affect 
their phase coherence and in order to implement them as building blocks for 
quantum computation schemes, it is necessary to understand which sources of 
decoherence affect the short time dynamics of flux qubits and reduce as much 
as possible their effect. 

Several sources of dissipation for flux qubits have been discussed through- 
out the literature [100], background charge fluctuations (r^ 0.1s), as well 
as quasiparticle tunneling in the superconductor with a non-vanishing subgap 
conductance (r^ ~ 1ms). Nuclear spins in the substrate have also been consid- 
ered as a possible source of dissipation. Static random magnetic fleld produced 
by the nuclear spins may induce shifts in qubit frequencies, but no dephasing 
is expected until a typical nuclear relaxation time, which can be very long, up 
to minutes, due to the slow dynamics of nuclear spins. 




(115) 



(116) 
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However, the most efficient source of dissipation for flux qubits is represented 
by fluctuations in the external circuit that produce fluctuating magnetic fluxes 
through the loop. The coupling of flux degrees of freedom of the qubit to the 
dissipative environmental elements is well described in the graph formalism 
described in [96]. In the Born-Markov approximation, it can be shown that the 
Redfield tensor, written in the eigenbasis {|n)} of TLs, is entirely determined 
by 



Rer!:L = • <^)/™(m . <^)nfc'^(l-n.|)-^^^^, (117) 

2 J(uj') 

ImrJ+^fc = -(m • " 'fUJil^nkD-P / dco ^ \ (118) 

X (u - UJnk coth , (119) 



where (3 = ksT and the m • if appears in the Hamiltonian for the system-bath 
coupling, m being related to the topology of the dissipative circuitry. 

In the two-level approximation, the rates for transitions from the two-level 
subspace to higher states can be neglected, and Eqs. (jll7p and (jllOp simplify 
and the dynamics of the 2x2 density matrix of the system can be cast in 
the form of Bloch equations for the dynamics of a pseudo spin 1/2. In this 
framework the relaxation matrix contains just two rates, T{~^ and for the 
decay of the longitudinal and transverse pseudo spin component, respectively. 
The latter in turn is limited by relaxation time Ti and pure dephasing time 
T^, I/T2 = I/2T1 + 1/T^, and the two rates are given by [96] 



;l = 4|(0|m.<^|l)|V(u;oi)coth-^, 
Ti 2kBT 

J(cu) 



;^ = |(0|m.(^|0)-(l|m.(^|l)| 



2kBT. 



(120) 
(121) 



Typically, T<^ can be made to diverge for an appropriate choice of external 
fluxes such that (0|m • (p\0) = (l|m • ip\l). However, this divergence is not ex- 
pected to show up experimentally, since it will be cut off by other mechanisms. 
The two lowest energy states, eigenstates of the Hamiltonian Eq. (I116p . are 
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given by 




(122) 



1) 



V2 



1 



f./l-— |L) - /l + — 

VV "^01 V ^01 / 



(123) 



where wqi = \/e^ + A^. Approximating the locahzed flux states |L) and as 
Gaussians centered in the minima of the double well, the relaxation rate Tj~^ 
and the pure dephasing rate TT^ are given by 



where S = {L\R) is the overlap between the two Gaussians. The vector A(/? 
connects the two minima of the double well. These relation are valid in the 
Markov limit and at in Born approximation, where the system-bath interaction 
is considered only at first order. By inspection of the previous formula it is 
clear that a symmetric double well potential, for which e = 0, let the dephasing 
time to diverge. This is realized for a value of the external applied magnetic 
flux $ext = ^0/2, being e oc ($ext/<l>o - 1/2). Moreover for m • A<^ = 
the environment is decoupled from the system, and both the relaxation and 
dephasing time diverge. 

In Ref. [96] an estimate of the leakage rate due to transition from the qubit 
states /c = 0, 1 to higher energy levels n =, 2, 3, . . . . outside the qubit subspace 
can be quantified from Eqs. ()117p and (|119p . 



In the regime 77 S> A, 5, e, where r/ is the energy splitting between the lowest 
two states | L) and | R) and the third energy level, and 5 is the coupling between 
the qubit subspace and the next higher level, the dominant leakage occurs with 
rate 




(124) 




(125) 




(126) 




(127) 
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and in the regime in which the two-level approximation is well defined, A 3> 
ksT, thermally activated leakage is strongly suppressed. 

1.20.2 Thermal photon noise induced dephasing. Besides magnetic flux 
fluctuations, an important source of dephasing is represented by thermal pho- 
ton noise in the read-out part of the circuit. To measure the state of the qubit a 
superconducting quantum interference device (SQUID) is coupled to the qubit 
via mutual inductance. When the SQUID is biased by a current that has a 
value smaller than the critical current, the SQUID acts just as an effective 
inductor, whose linear inductance depends on the qubit current state. The 
critical current at which the SQUID switch to a normal state can have two 
values /i*^^ and , according to the two qubit possible current states. A bias 
current pulse of amplitude Ib , < Ib < Ic''^ , allows to discriminate the two 
qubit states. The read-out apparatus consisting of a dc-SQUID and a shunt 
capacitor Cg form a weakly damped harmonic oscillator of frequency (Jhoj that 
is detuned from the qubit frequency uoq. The presence of n photons in the 
harmonic oscillator induces a shift in the qubit frequency, ujq^n — ^q,0 = nSuQ, 
where the shift per photon depends on the effective qubit-oscillator coupling. 
Assuming that the pure dephasing time is much larger that the inverse of 
the damping rate k, ^ thermally excited photons in the oscillator 
produce a dephasing [86, 101, 102] 



where n = {cx:p{hiJ\io / UbT) — 1)^^ is the thermal average number of photon in 
the oscillator. A similar effect has been observed in an experimental work [87] 
in which a charge qubit is coupled to a superconducting waveguide resonator, 
slightly detuned from the qubit frequency. There, opposite to the case here 
described, the oscillator is driven and a shift and a broadening in the qubit 
resonance frequency appears, as a consequence of an ac-Stark shift and of 
photon shot noise. 

The flux qubit of Ref. [102] has been engineered with four Josephson junction 
to ensure a symmetric qubit-SQUID coupling [103]. In the usual design of the 
Delft qubit [92], the two symmetric arms of the SQUID render the qubit 
immune to bias current Is fluctuations. At zero dc bias, = 0, a small 
fluctuating current caused by the finite impedance of the external controls 
is divided equally into two branches of the SQUID loop and no net current 
flows through the three-Josephson junctions of the qubit line. However, the 
double layer structure of the Josephson junctions, being an artefact of the 
shadow evaporation technique used to construct Josephson junctions, induces 



(128) 



r0 = 



n(n + l)(5a;j 



'o 
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asymmetry in the circuit. Using a forth much larger Josephson junction, for 
which the Josephson energy can be usually neglected, symmetry in the double 
layer structure is restore and effects of fluctuations in the bias current are 
suppressed. 

The qubit energy bias can be written as the sum of two contributions, e = r]+ 
A, where rj = 2/p($ext — $o/2) {Ip is the qubit persistent current) is controlled 
by the external flux $ext and A = 2IpMJ{lB)/h which depends on Ib via the 
SQUID circulating current. This dependence has two crucial consequences: 
first the qubit bias point $*xt which dug/d^ext = results shifted by the 
measurement pulse. Therefore it is possible to operate the qubit at the flux- 
insensitive point, while keeping a difference in the expectation value of the 
current in the two qubit states, which is a crucial requirement for measuring the 
qubit state. Second a coupling between the qubit and the external harmonic 
oscillator, the so-called plasmon mode, arise, with an interaction Hamiltonian 

Hq-ho oc [gi{lB){a + at) + g2{lB){a + a^fja,, (129) 

where gi{lB) oc (dX/dlB) and 92(13) oc (cfX/dl'^) [101]. For a particular 
that realizes dX/dlB = 0, it is possible to switch 51 (/s) off [102,103]. 

Working at these optimal point, the c|ubit is immune from external 
flux and bias current fluctuations at first order, dtOq/d^exti^l^t^ -^b) ~ 
du>q/dlB{^lxti ^b) — 0- '^^^ ^^ift P®^ photon of Suq is given, at second or- 
der perturbation theory in Hq-ho [101,102], by 



Sloo = 4 



{gi (Ib) sin df ^ ^'^ - g2 {Ib) cos 6 



(130) 



where cos 6 = e/uiq. For some value e*{lB) < one obtains Sloq = 0. In 
Ref. [102] , via spectroscopy the authors demonstrated the existence of a line 
€*{Ib), that includes Ib = Ib and e = 0, providing an optimal point with 
respect to bias current noise, flux noise, and photon noise. Measurements of 
the qubit spectral line shape at the optimal point showed, for the particular 
sample Ref. [102], a twin peak structure, which could arise from strong coupling 
to one microscopic fluctuator. An effective dephasing time t2 = 2/tt{wi + W2) 
is obtained by fitting the peaks with two Lorentzians of width wi and W2- 
Measurements of the spin-echo decay time Tcchoi particularly indicated in case 
of relatively high frequency noise, as photon noise in the plasma mode that 
occurs at K 130MHz, gives at the optimal point Techo = 3.9//S. By studying 
the variation of Techo and t2 as a function of e, a sharp peak is found at e = 
for Ib = Ibi while for Ib = the peak shifts towards e < 0. The variation of 
the maximum in t2 as a function of Ib show that the maximal coherence time 
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Figure 9. a) Schematic representation of the phase qubit circuit, constituted by a current-biased 
Josephson junction, b) Anharmonic potential U{5) showing the two lowest energy states |0) and 
|1), separated by an energy splitting Eqi- The plasma frequency ujp is given by the local quadratic 
curvature of the potential at the bottom of the well, and Af7 is the potential barrier that separates 
the two energy levels in the well from a continuum of levels on the right side of the barrier. 



is not obtained at e = 0, as it would be expected for flux or bias current noise. 
On the other hand, it fits with the line (*{Ib) for which 5ujq = 0, suggesting 
that thermally induced photon noise, rather than flux noise or bias current 
noise, is responsible for the qubit dephasing. For a temperature T = 70mK 
and quality factor Q = 150, which yields a mean photon number n = 0.15, the 
dephasing time Eq. (|128p closely matches the spin-echo measurements. 



1.21 The superconducting phase qubit 

The phase qubit works in a regime in which Ej ^ Eq and the circuit consists 
of a loop with a single large Josephson junction, as shown in Fig. [9] (a). The 
circuit is biased with a current / typically driven close to the critical current 
/q of the junction. The Hamiltonian of the system is 

where ^0/21: = h/2e is the superconducting flux quantum. Charge and phase 
operators, Q and (5, that correspond to the charge on the Josephson junc- 
tion capacitance and the superconducting phase across the Josephson junction 
respectively, are conjugate variables that satisfy the canonical commutation 
rule [5, Q] = 2ei. For a large area junction Iq(^q/2it = Ej ^ Ec = e^/2C 
the superconducting phase has a well defined value and quantum mechanical 
behavior can be observed. The Josephson inductance and the junction capac- 
itance form an anharmonic "LC" resonator. The anharmonic potential as a 
function of the superconducting phase across the junction can be approxi- 
mated by a cubic potential parametrized by the potential barrier AU{I) = 
(2y^2/o^o/37r)[l — I/Iq]'^^'^ and a classical plasma oscillation frequency at the 
bottom of the potential well ujp{I) = 2i/^(27r/o/^>oC)^/2[;L _ ///o]i/4_ The two 
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qubit states |0) and |1) are encoded in the two lowest quantum states in the 
potential well, and have energy splitting Eqi = hjjjp{l — bhup/SGAU). Unlike 
in a flux or charge qubit, in a phase qubit the state (|0) or |1)) is exclusively 
distinguished by the phase wavefunction, and not by any macroscopic quan- 
tity, such as current or charge. Transitions are driven by applying microwaves 
at frequency ojqi = Eqi/K. For more details on the phase qubit we refer to 
Ref. [94]. 

Coherent control of the qubit is obtained through the bias current 

I{t) = Idc + h/fit) + IfMwcit) cos(u;oit) + I^iwsit) sin(a;oit), (132) 

where Ii/f, Ifiwc, and Ifj,ws are varied in time slowly compared to 2tt/ujqi. In 
the frame rotating with frequency luqi, the qubit Haniiltonian is given by 

where e = y^^^T^i) crx,y,z are the Pauli operators. 

1.21.1 Decoherence of a phase qubit due to an arbitrary noise source. 

Since the qubit is controlled by the bias current, noise in the bias current can 
represent a source of decoherence for the qubit. In Ref. [104] a physical picture 
of decoherence is presented for a phase qubit. 

In a Bloch picture, the state of the qubit is represented by \^) = 
cos(^?/2)|0) + e**^ sin(^/2)|l). Low frequency noise induces fluctuations is the 
longitudinal z component of the pseudo-spin representing the qubit that lead 
to dephasing of the qubit. The phase noise after time t is 

Ut) = ^ fdt'lnit'), (134) 
C^dc Jo 

and it arises from current noise In{t)- The magnitude of the phase noise is 
described by (^^^(i)), and it can be obtained through the noise spectral density 

(fir. \ 2 /'uioi/2n 
dfSi{f)Wo{f), (135) 

where Sj{f) is given by the Fourier transform of noise correlator {Init)In{0)) = 
dfSi{f)cos{2TTft), the spectral weight Wo{f) = sin^i-K ft) / {-k ff , and the 
integral on the frequency has been cutoff for frequency greater than t(;oi/27r. 
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This last assumption is justified since for tliose frequencies the noise current 
flows mainly through the junction capacitance, rather than the junction it- 
self, thus not substantially affecting wqi- The magnitude of noise is defined 
as the mean-square amplitude of the current noise at frequency / per 1 Hz 
bandwidth. For low frequency / < 1/t, VFo(/) is rather constant, whereas it 
decreases as 1//^ at higher frequencies. As a consequence, phase noise affects 
the qubit dynamics only at low frequencies for most noise sources. For constant 
(white) noise 5^, one has 

At higher frequencies close to wqi, noise induces transitions between the two 
qubit states |0) and The current that controls these transitions is given by 
I nwcif) ^^^{^oit) + -^fi«is(i) sin(u;oi^)) and mixing from noise around frequency 
ujQi can be understood as low frequency noise in I^wc{t) and Ifj_ws{t)- Random 
fluctuations along the transverse x and y components of the qubit induce 
transitions between the qubit states. For constant spectral density around wqi 
given by 2S'/(a;oi/27r), an application of the previous results gives 

{Ol{t)) = (0^(0) = ^ 5,(a;oi/27r)t. (137) 

The random angles x = 4>,()x and 9y are assumed to be Gaussian dis- 
tributed, with zero mean and mean squared noise (x^) previously calculated, 
dp{x)/dx = exp(-xV2(x'))/V'2^. 

If the qubit is initially in the ground state, that is when it is parallel to the 
z direction in the Bloch sphere, when the noise is small, it is immune to phase 
noise at low frequency. However, transverse noise around frequency wqi can 
induce transitions in between the qubit states. The probability po to be in the 
state |0) is given in the Bloch picture by cos^(0/2) ~ q,os^[.^6% + 6^2]. With 

the values previously obtained for the mean-square noise, averaging over the 
Gaussian distribution gives 

The rate 1/Ti describes absorption and emission rate for the stimulated tran- 
sitions 0—^1 and 1 — > 0. Since low frequency noise cannot add energy Hloqi, 
there is no contribution from phase noise. 
Effects of noise on a superposition state can be understood within a "Ram- 
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sey fringe" picture. Through SlTt/2 pulse the qubit is rotated from the ground 
state |0) to the state (|0) + |l))/\/2, that points in the x direction in the Bloch 
sphere, and left evolving for a time after which a subsequent 7r/2 pulse is 
performed and the qubit state is finally measured. During the evolution be- 
tween the two 7r/2 pulses, the state of the qubit can change due to noise in cf) 
and 9y, therefore, both phase and stimulated transitions noise affect the qubit 
dynamics. In this case the total decoherence rate is given by the Korringa rela- 
tion [9] I/T2 = l/T^ + l/2ri, where l/T^ is directly extracted from Eq. (fT36]l 
l/T^ = idu;oi/dh)^Syi. 



1.22 Decoherence due to 1/ f noise 

Much effort has been spent recently to understand how noise at low frequen- 
cies affects the dynamics of superconducting qubit, both from a theoretical 
and an experimental point of view. In particular, signatures of non-Markovian 
dynamics are believed to be due to 1 // noise. Both charge and flux 1 // noise 
contribute to decoherence, the former affecting mostly the dynamics of su- 
perconducting qubits based on the charge degrees of freedom, and the latter 
affecting flux and phase qubits. Here we concentrate only on 1// flux noise, 
and refer to the literature for the case of 1/f charge noise [105,106]. 



1.22.1 Model for 1/f flux noise in SQUIDs. The origin of 1// flux noise 
in low- Tc devices has not yet been completely understood. Critical-current 
fluctuations in Josephson junctions are believed to arise from trapping and 
release of electrons in traps in the tunnel barriers. In Ref. [107], a model for 
1/f noise in low- Tc devices is proposed, based on the assumptions that it arises 
from thermally activated hopping of unpaired electrons on and off defects. The 
spin of the electron is locked in direction while the electron is trapped, and 
the directions are randomly change from trap to trap. Uncorrelated changes 
of these spin directions give rise to random telegraph signals that produce a 
1/f power spectrum. The electron is assumed to occupy the low-energy spin 
direction during the entire time it resides on the defect. For zero magnetic 
field B, transitions between the two degenerate Kramers' doublet are strongly 
suppressed, implying that direct phonon scattering is forbidden. On the other 
hand, the magnetic field is not strictly zero, and fluctuating dipole fields can 
arise from neighboring defects. 

The magnetic moment of a defect M = ^^(L + 2S) can be locked as a 
consequence of spin-orbit coupling, that makes it stable with respect to these 
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weak fluctuations. The locking effect can be modeled by the Hamiltonian 

H= Yl +AL-S + MsB- (L + 2S). (139) 

i=x,y,z 

The unpaired electron occupies a p orbital, and Vx^y^z are the matrix elements 
of the crystal-field potential. The spin-orbit coupling constant, depending on 
the different kind of defects, can vary in magnitude within a large range of 
values, but for defects involving atomic weights near that of Si, |A| « 300 K. 
The crystal-field Vx^y^z can be at most ^ 2000 K. In defects for which A < 0, 
L and S are parallel and M is large, while for defects for which A > L and 
S are antiparallel, and M is close to zero. Therefore, the A < defects are 
expected to be most important for flux noise. 

A random distribution of defects over the substrate is assumed, and the flux 
noise coupled into a SQUID by a spatially random distribution of electron 
spins, whose orientations fluctuate, is calculated. The magnetic moment is 
represented by a small current loop that couples to the SQUID loop by a 
mutual inductance M{x,y). The SQUID loop is schematized like a square 
frame of inner and outer dimensions of 2d and 2D, and thickness W = D — d, 
lying in the plane z = 1 /xm. The current loop can lie in the z = /zm plane, 
resulting in a perpendicular moment, or in the x or y plane centered at z = 0, 
resulting in a in-plane moment. The small loop current has an effective area 
A = (0.1 /[zm)^, and a current i flows in it, such that Ai = hb, with flux 
per Bohr magneton given by ^s/fJ-B = M{x, y)/A. Perpendicular and in-plane 
flux per Bohr magneton show a qualitative opposite behavior as function of 
the position in the SQUID loop plane, the former peaking on the edges of the 
superconductor, and the latter peaking at the midpoint of the superconductor, 
both falling off away from these points. The total mean square normalized flux 
noise coupled to the SQUID is given by 

pL+D f-x m'^(t 

ms?) = Snfil dx dy ^' (140) 

where n is the areal density of defects, and the integral over the plane is cutoff 
a distance L away from the SQUID loop. The mean square noise is given by 

{{5^st?) = m^si,x?) + msi,y?) + ((5$sp)')]/3, and the spectral density 
S'$(/) = a/f is extracted by {{5^stf) = ajj^df/f = aln(/2//i), and, for 
/2//1 ~ 10^3^ 

S^if)/^l ~ ms/^of)/m- (141) 
Noise levels in agreement with the observations are obtained for n ~ 
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10^^ m ^. However, they strongly vary with the geometry of the the SQUID 
and the tunnel barriers, and with the fabrication details. 

1.22.2 Decoherence of flux qubits due tol/ f noise. Recently experiments 
[108,109] have reported the behavior of the echo signal in flux qubits at various 
bias conditions. The energy splitting depends on the applied external magnetic 
flux, and thus it is sensitive to flux noise. As a result they found that the at 
the optimal point, where the energy splitting is insensitive at first order to 
magnetic flux fluctuations, the coherence time is limited by energy relaxation 
processes, Ti, and the dephasing of the flux qubit is mostly determined by the 
high-frequency noise S^{uj ~ A/h). 

In a Markovian approximation the energy relaxation contributes to the de- 
phasing process via l/TI'^'^" = l/2Ti + F^, where is the pure dephasing 
rate. For dephasing dominated by magnetic flux noise with a smooth spectrum 
near uj = 0, the pure dephasing rate is given by F* ~ {dAE/d^cxt)^S{uj = 0), 
where AE = ye^^^^^^J^TpA^. Taking into account the fact that, from their ex- 
perimental data, the relaxation time Ti is almost independent on the external 
applied magnetic flux, and that close to the optimal point the dAE/d^cxt 
^ext) a parabolic behavior of 1/T2'^^° is expected in this region. On the other 
hand, away from the optimal point, they observed a linear increasing of l/T2'^^° 
with respect of the applied external magnetic flux. Therefore, the experimental 
data cannot be explained within the framework of Bloch-Redfield decoherence 
theory, in which the assumption of short time correlated noise (white noise 
around a; = 0) holds. 

The experimental observations can be explained with the presence of 1// 
flux noise. In Refs. [108,109], in order to separate the contribution to dephasing 
due to direct transitions between energy levels, the echo signal is expressed 
as p{t) = e~*/^-^^/9echo(i)) with p{t) = (cr^). In this case the expectation value 
of (T^ is given by a non-exponential decay curve, and at the end of the echo 
sequence 



where = {^o/fi)dAE/d^ext- Assuming Gaussian statistics of the fluctua- 
tions of the external flux, the decoherence rate can be expressed through the 
noise spectral function, S'$(a;) = (I/tt) dt cos ut{^{t)^{0)), 




(142) 



p{t) = exp - 



12 poo 



sin^ {cot /A)' 
(a;i/4)2 



(143) 
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For a 1// spectrum = / uj one obtains 



Pecho(i) = exp [-{T^tf] 



Tl = \v^WA^\n2. 



(144) 



In Ref. [110], motivated by the experimental results of Refs. [108, 109], a 
theoretical analysis of the dephasing due to 1// noise in flux qubits has been 
presented. The problem is described by choosing spin-fluctuators as source 
of low-frequency noise. Random switching of a fluctuator between its two 
metastable states gives rise to a random telegraph noise. Transitions in fluctua- 
tors with energy splitting larger than the temperature are strongly suppressed 
and only fluctuators whose energy splitting is smaller than the temperature 
contribute to the qubit dephasing. For a random telegraph process, ^(t), in 
which a switching between the values ±1/2 takes place at random times, the 
probability to make n transition in a certain amount of time r is Poisson dis- 
tributed. As a consequence, the time correlation function ('^(i)^(O)) = e~^^^^ /A 
decays exponentially, with I7I given by the rate of the transition between the 
to states of the fluctuator. The contribution to the noise spectrum of a random 
telegraph process, given by the Fourier transform of the correlation function, 
is a Lorentzian, ^ /A'k{uJ^ + 7^). The effects of many uncorrelated fluctuators 
coupled to the qubit via constants Vi simply add. Considering a large number 
of effective spin-fluctuators and assuming no correlations between couplings Vi 
and switching rates 7^, the noise spectral function is give by 



The distribution ^^(7) depends on the details of the interaction between the 
fluctuators and the qubit. Following [110,111], 7^7(7) ~ (71/7)0(7 — 70), where 
7o represent the maximal relaxation rate for fluctuators with a given energy 
splitting E, and A gives the amplitude of the 1// noise. It follows for the noise 
spectral function 



The spin-fluctuator model reproduces 1// noise power spectrum for u; ^ 70, 
but it predicts a crossover from a uj~^ to uj~'^ behavior, consequence of the 
assumption of a maximal switching rate 70. Expressing the fluctuation of the 
magnetic flux as sum of the contributions of the independent fluctuators <I>(t) = 
Yli ill a Gaussian approximation, substitution of Eq. (|146j) in Eq. (I143|) 




(145) 




(146) 
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11.S9 11.60 11.61 11.62 

Bias Current I [jxA] 

Figure 10. Measured probability of state "1" versus microwave excitation frequency u}/2tt and bias 
current / for a fixed microwave power for a phase qubit, obtained from Ref. [112]. Dotted vertical 
line indicate spurious resonators. (With permission from the Authors). 



leads to 

r,^-lnpHWx{^f/. (147) 

Therefore the echo signal is expected to decay in a Gaussian way at long times 
t » 7(^^, while at short times t <^ '^q'^ a faster decay is expected, Kg oc t^. 

Using the method of stochastic differential equations (see [110] and refer- 
ences therein), an estimate for the non-Gaussian case provides for the loga- 
rithm of the echo signal /Cgf [110] 

r -i^AvH'^/Q , t « 7o-i 

/Crf ^ I ln2AvH^ , « t < (148) 
[ avAt , v^^ <C t, 

where v is the center of a sharp peaked distribution of the couplings between 
the fluctuators and the qubit, and a « 6. A new decaying behavior arise in 
the case of large time t ^ v~^, that drastically differs from the predictions 
of Gaussian statistically distributed magnetic flux fluctuations. The reason 
laying in the fact at short time only fast "fast" fluctuators contribute to the 
dephasing, giving rise, for v <^ j, to Gaussian decay, while at long time non- 
Gaussian behavior appears. 
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1.22.3 1// noise in superconducting phase qubit. In this section we re- 
view the latest results on decoherence of superconducting phase qubits [93,94]. 
Recent experiments [112,113] have pointed out that a dominant source of de- 
coherence for the phase qubit is represented by low frequency 1// noise, that 
is believed to arise from two level systems (TLS) in the insulating barrier of 
the tunnel junction as well as in the dielectric material surrounding the cir- 
cuit. In Fig. (Tola measurement of the transition frequency cjqi between the two 
qubit state, as a function of the bias current / and the microwave excitation 
frequency iv/2TT, shows a qubit line in which a number of spurious resonators 
appear, characteristic of energy-level repulsion predicted for coupled two-state 
systems. Near the resonators, the Rabi oscillations show beating, loss and re- 
covery of the oscillations with time, and rapid decrease of coherence ampli- 
tude. The beating behavior is consistent with the interaction of a qubit with 
another two-level system, but not with harmonic oscillator modes in the read- 
out SQUID. Moreover, each qubit has its own set of resonator frequencies and 
strength, indicating that the TLS have a microscopic origin. 

A new method to measure 1/f noise in Josephson junction qubits has been 
recently presented in [114] . It uses the resonant response of the qubit to directly 
measure the spectrum of the low-frequency noise, and allows to distinguish 
between flux and critical-current fluctuations by comparison of the noise taken 
at positive and negative bias. Remarkably it can yield low- frequency spectra 
below 1 Hz. Dephasing is produced by low-frequency fluctuations in the qubit 
energy, which in this study are believed to arise from magnetic flux noise in 
the qubit loop, with a spectral density that scales inversely with frequency, 
1/f. It turns out that the is flux-like noise predominates over critical-current 
noise. 

The possibility that flux noise is due to TLS defects in the native oxides 
of the superconductive film, as proposed in [107], is examined in Ref. [114]. 
Following [115] for a standard TLS model [116] a theoretical estimation of the 
flux noise spectral density, for a realistic geometry of the circuit loop, gives 
^^(l Hz) 10-3(^$o)2/Hz, about 4 orders of magnitude smaller than the 
measured flux noise. This estimate is based on the assumption that TLS fluc- 
tuations randomizes the defect magnetic moment; assumption highly question- 
able because TLS defects in typical oxides are not considered to be magnetic. 
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